V- 


u 


DOCUMINT  CONTROL  DATA  .RAO 

*ktty*ti¥  t  Itta/tf  ttton  of  fill#,  b*dy  •/  ilflMir  0*4  inditing  m nitttlm  matt  bt  mnfwrtd  whtn  th$  artrttl  ttmff  ft  rfottfUidl 


I  O^iOlNAfiNO  «C  fivifv  (C  try  art  It  tulhtf) 

Princeton  Univorsity 

Department  of  Aerospace  &  Mechanical  Science! 
Princeton,  New  Jersey  08540 


I  RKRORT  T  l  T  i.  4 

Micro-Rocket  Impulsive  Thrusters 


u.m^out  fccunirv  et 

Unclassified 


'*Ten'ti*«o  pffii" 


UTMOAifi  (Flm  ntm,,  ml 


tpo'l  «ntf  Intluilvi  It  ft) 


til  ntmt) 


Leonard  H.  Caveny  and  Martin  Summerfield 


•  RIRORT  OA  Tf 

November  1971 


It,  C  ON  TRACI  ON  OWANT  NO 

DAHC  60-71-0 0034 

b.  *mojic  t  no 


70,  TOTAL  NO.  OF  PAGII  7b.  NO.  OF  RIFS 

100  49 


90,  OftfOINATO***  RCRORT  NUMIM(SI 

AMS  Report  No.  1014 


OTNIR  RKROR T  NOIH  (Any  othtr  numbtrt  Uttf  mt/  bt  tttlgntd 
thlt  rtporl) 


\0  Of  •  T  Rl  SU  T  ION  STATEMENT 

Distribution  of  this  document  is  unlimited 


It.  IPONIORINQ  MILITARY  ACTIVITY 


U.  S.  Army  Advanced  Ballistic 
Missile  Defense  Agency 


Light  weight  solid  propellant  motors  which  supply  high  thrusts 
(>500  Ibf)  for  durations  on  the  order  0.010  sec.  (referred  to  as 
impulsive  thrusters)  require  special  analysis  of  items  such  ac 
internal  ballistics,  propellant  combustion,  ignition  stimuli,  exhaust 
plume  envelope,  and  inert  hardware.  An  almost  explosive  ignition 
and  pressuri cation  coupled  with  the  requirements  for  reproducible 
thrust  versus  time  programs  are  analyzed  by  a  mathematical  model  that 
emphasizes  gas  flow  and  propellant  combustion  dynamics.  Studies  were 
directed  at  analyzing  a  center-vented  motor  approximately  7  inches 
long  and  0.3  inches  in  diameter  with  a  total  impulse  of  3.6  lbf-sec 
and  a  web  time  of  0.006  sec.  A  survey  of  the  rocket  motor  concepts 
revealed  that  the  desired  thrust  versus  time  program  can  be  achieved 
by  internal  burning  grains  and  existing  high  burning  rate  composite 
propellants.  Limitations  were  placed  on  volumetric  loading  density 
as  a  result  of  uncertainties  of  how  the  center-vented  flow  affects 
thrust  vector  alignment  and  nozzle  flow  losses.  The  mathematical 
model  for  transient  internal  ballistics  developed  during  the  study 
can  be  applied  to  a  wide  range  of  high  performance  rocket  motors 
which  experience  rapid  ignition  and  pressurization  transients. 


DD  .’.r.,1473 


Unclassified 


Best 

Available 

Copy 


Unclassified 


*-101  - 


a*tr  I  of  i  e»*k,a 


Solid  propellant 
High  pressure  rocket  motors 
Propellant  combustion 
Internal  ballistics 
Rocket  motor  ignition 
Very  small  rocket  motors 


Unclassified 


Security  Clc»niflc»*ie* 


MICRO- ROCKET  IMPULSIVE  THRUSTERS 
by 

Leonard  H.  Caveny  and  Martin  Summerfield 
AMS  Report  No.  1014 

November  1971 


Performed  under  Contract  No.  DAKC  60-71-C-0034  issued  by  the 
U.  S.  Army  Advanced  Ballistic  Missile  Defense  Agency. 


Qualified  requestors  may  obtain  additional  copies  from  the 
Defense  Documentation  Center. 


Guggenheim  Laboratories  for  the  Aerospace  Propulsion  Sciences 
Department  of  Aerospace  and  Mechanical  Sciences 


PRINCETON  UNIVERSITY 
Princeton,  New  Jersey 


j — DISTKIBOTiON  statement  a_ 
for  public  rel's-se; 
DistriV’.ik-n  Unjr.’'ifcd _ 


FORWARD 


This  is  the  final  report  issued  under  Contract  No.  DAHC 
60-71-C-0034.  Reported  herein  is  work  performed  over  an 
eleven  month  period  starting  15  November  1970,  by  the 
Aerospace  and  Mechanical  Sciences  Department  of  Princeton 
University. 

The  work  was  monitored  Dr.  David  C.  Sayles  of  the 
Huntsville  Office  of  the  U.  S.  Array  Advanced  Ballistic 
Missile  Defense  Agency. 


ABSTRACT 


Light  weight  solid  propellant  motors  which  supply  high 
thrusts  (>500  Ibf)  for  durations  on  the  order  0.010  sec. 
(referred  to  as  impulsive  thrusters)  require  special  analysis 
of  items  such  as  internal  ballistics,  propellant  combustion, 
ignition  stimuli,  exhaust  plume  envelope,  and  inert  hardware. 

An  almost  explosive  ignition  and  pressurization  coupled  with 
the  requirements  for  reproducible  thrust  versus  time  programs 
are  analyzed  by  a  mathematical  model  that  emphasizes  gas  flow 
and  propellant  combustion  dynamics.  Studies  were  directed 
at  analyzing  a  center-vented  motor  approximately  7  inches 
long  and  0.3  inches  in  diameter  with  a  total  impulse  of  3,6 
lbf-eec  and  a  web  time  of  0.006  sec.  A  survey  of  the  rocxet 
motor  concepts  revealed  that  the  desired  thrust  versus  time 
program  can  be  achieved  by  internal  burning  grains  and  existing 
high  burning  rate  composite  propellants.  Limitations  were 
placed  on  volumetric  loading  density  as  a  result  of  uncertain¬ 
ties  of  how  the  center-vented  flow  affects  thrust  vector 
alignment  and  nozzle  flow  losses.  The  mathematical  model 
for  transient  internal  ballistics  developed  during  the  study 
can  be  applied  to  a  wide  range  of  high  performance  rocket 
motors  which  experience  rapid  ignition  and  pressurization 
transients . 
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INTRODUCTION 


In  some  solid  propellant  applications  it  is  desirable  to 
produce  a  controlled  high  thrust  for  a  very  short  time,  that 
is,  to  make  the  burn-out  time  no  longer  than  a  few  milli¬ 
seconds.  All  of  this  must  be  accomplished  in  high  perfor¬ 
mance  systems,  i.e.,  maximum  total  impulse  per  motor  weight. 

We  will  refer  to  such  systems  as  impulsive  thrusters  or  micro¬ 
rockets.  Typically,  these  systems  use  less  than  one-tenth 
pound  of  propellant  and  are  less  than  one  inch  in  diameter. 
Broadly  speaking,  two  internal  ballistics  approaches  are  pos¬ 
sible.  The  first  approach  is  to  desian  the  propellant  grain 
so  that  the  effective  web  thickness  of  the  propellant  grain 
is  very  small  and  to  take  advantage  of  increased  burning  rates 
at  high  combustion  pressures  (circa  15,000  psi) .  However, 
very  thin  webs  are  practical  only  if  the  unsupported  web  thick¬ 
ness,  r  ,,  is  an  appreciable  fraction  of  the  propellant 
diametereB  (e.g.,  2Twey/DDror)  >  0.1)  and  acceleration 

and  gas  dragforces  are  small?  °Because  of  the  difficult  prob¬ 
lems  in  grain  design  in  pursuing  this  route,  a  second  approach 
may  be  more  attractive.  The  second  approach  is  to  utilize 
high  burning  rate  propellants,  which  may  be  faster  by  five  or 
more  times  than  the  standard  propellants.  An  additional 
consideration  to  the  development  of  impulsive  type  microrockets 
is  establishing  the  most  efficient  case,  nozzle  and  igniter 
designs. 


This  report  summarizes  the  investigations  (conducted  at 
Princeton  University)  to  determine  and  evaluate  approaches  to 
the  development  of  solid  propellant,  impulsive  type  micro¬ 
rockets.  The  applications  of  these  investigations  were  di¬ 
rected  at  satisfying  the  Hit  propulsion  requirements.  The 
study  at  Princeton  University  consisted  of  the  following  items: 

I.  Analysis  of  approaches  to  achieve  controlled  thrust 
for  very  short  times  and  the  development  of  a  com¬ 
prehensive  mathematical  model  for  predicting  the 
performance  of  micro-rocket,  impulsive  thrusters. 

II.  Analysis  of  methods  of  achieving  ultra-hiqh  burn¬ 
ing  rate  propellants. 

III.  Study  of  operating  characteristics  and  methods  of 
maximizing  performance  of  the  most  promising  im¬ 
pulsive  thruster  concepts. 


The  major  objective  of  this  study  is  to  establish  a  sound 
technical  basis  for  evaluating  the  various  competitive 
approaches  to  the  impulsive  thrusters  for  Hit. 


SUMMARY  OF  REQUIREMENTS 


Before  beginning  the  discussions  of  the  approaches  to  im¬ 
pulsive  thrusters,  a  discussion  of  the  operating  conditions 
will  help  to  focus  more  clearly  on  the  requirements  that  a 
practical  motor  design  and  propellant  must  meet.  Rapid  ignition 
and  pressurization  (all  within  a  few  milliseconds  or  less)  are 
characteristic  requirements.  Thus,  the  propellant  should  be 
capable  of  withstanding  an  almost  explosive  ignition  impulse. 

As  in  most  solid  rocket  motors,  it  will  be  necessary  to  pro¬ 
vide  pressure  level  control  by  means  of  a  prescribed  burning 
surface  area  versus  web  thickness:  thus  the  high  burning  rate 
propellants  should  lend  themselves  to  being  molded  or  cast 
into  various  geometries.  Unfortunately,  some  of  the  mechanical 
approaches  to  achieving  high  burning  rates  are  better  s'uited 
for  end-burning  grains  rather  chan  internal -burning  grains. 

Of  major  importance  is  that  the  process  which  is  used  to  in¬ 
crease  burning  rate  must  be  reproducible  and  repeatable  in 
practical  rocket  motor  configurations. 

Table  1  shows  in  summary  form  the  basic  requirements  of 
the  Kit  propulsion  system  (a  complete  specification  is  given 
in  Table  3) .  a  pressure  versus  time  program  corresponding  to 
the  Hit  requirements  is  shown  in  Fig.  1.  Note  the  very  short 
burning  time  and  the  high  pressure  cf  the  system.  The  mass 
fraction  requirement  becomes  very  restrictive  when  it  is  re¬ 
viewed  in  terms  of  the  center-vented  requirement  and  miniature 
size  of  the  components.  The  photograph  of  a  model  of  the  Hit 
impulsive  thruster  shown  in  Fig.  2  gives  an  impression  of 
the  size  and  general  configuration-  Note  the  cross-section 
of  the  propellant  grain  showing  the  very  small  star  points 
and  port  passage.  Small  errors  in  aligning  the  internal  port 
of  grain  can  greatly  alter  the  burn-out  characteristics  of 
the  motor. 


TABLE  1 


SUMMARY  OF  REQUIREMENTS  FOR  HIT  SYSTEM* 


IT,web'  lbf  sec 

3.6 

fcweb'  sec 

<0.008 

Mass  fraction,  Wprop/WMQ 

Minimum  length,  L^,  in 

>0.49 

6.0 

Safety  factor 

1.4 

Peak  thrust  F  ,  lbf 

ITlo  X 

700 

Ignition  delay  to  first  thrust,  sec 

<0.0015 

Variation  of  centroid  of  pulse, 

3c,  sec 

<0.00125 

Chamber  is  center  vented  {i-e.,  side  vented) 

Thruster  is  rotated  on  the  outer 
circumference  of  a  cylinder  and 
its  nozzle  is  aligned  perpendicular 
to  the  cylinder  centerline. 

_ 

An  example  specification  is  gi^en  in  Table  3 


PRESSURE  AT  WHICH 


Pig.  1  Events  that  occur  during  the  operation  of  an 
impulsive  thruster  (representative  of  Hit 
type  of  application). 


configuration  (In  the  foreground  is  a  cross-section  through 
propellant  grain). 


SURVEY  OF  APPROACHES 

Over  the  last  two  decades,  many  attempts  to  develop  moder¬ 
ate  size  propulsion  systems  (between  2  and  4  inches  in  diameter) 
with  high  mass  discharge  rates  have  been  reported.  The  large 
part  of  this  work  has  been  directed  ac  increasing  the  effective 
burning  rate  of  the  propellant.  To  a  lesser  extent,  the  feasi¬ 
bility  of  various  methods  of  employing  very  thin  webs  has  been 
examined.  Because  of  the  requirement  for  high  total  impulse 
per  motor  weight  (IT/WM0)  systems  such  as  the  cap  pistol  motor 
were  not  considered.  The  high  thrust  requirement  eliminates 
from  consideration  systems  such  as  the  subliming  solids  and 
monopropellant  vernier  control  motors. 

The  various  methods  of  extending  burning  rates  include 
the  use  of  mechanical  modifiers!^  (metal  fibers  and  tubes), 
very  small  oxidizer  sizes,  burninq  rate  catalysts^,-  and  specially 
formulated  propellant  ingredients2* .  These  methods  have  not 
achieved  the  full  range  of  increases  in  buring  rate  that  are 
being  sought.  However,  some  of  the  requirements  for  moderately 
high  mass  discharge  systems  can  be  met  by  high  rate  propellants 
in  conjunction  with  high  surface  area  (i.e.,  thin  propellant 
webs)  reasonably  structurally  sound  propellant  grains.  As  the 
desired  operating  time  decreases  and  the  motor  case  diameter 
increases, there  are  limitations  in  following  the  approach  of 
using  high  burning  rate  propellants  in  conjunction  with  very 
thin  web,  internal  burning  propellant  configurations,  (e.g.,. 
a  wagon  wheel  type  of  configuration) .  Configurations  with 
small  Tweb/Dprop  ratios  will  not  withstand  the  loads  produced 
by  high  velocity  gas  flow  and  acceleration  of  the  system. 

When  the  fragile  propellant  star  points  are  broken  off  and 
ejected  through  the  nozzle,  erratic  pressure  traces  and  in¬ 
efficient  motor  operation  will  result.  Figure  3  shows  the 
extent  to  which  decreasing  operating  time  affects  the  thick¬ 
ness  of  unsupported  portions  of  propellant  configurations. 

For  the  particular  burning  rate  (10  in/sec)  which  was  con¬ 
sidered  in  Fig.  3,  reasonable  web  fractions  are  obtainable 
for  the  H'.t  operating  time  and  case  diameter.  However,  such 
high  burning  rates  are  available  only  in  recently  developed 
propellants  many  of  which  can  not  be  considered  state-of-the- 
art  with  respect  to  the  Hit  requirements. 

Variov  :  mechanical  methods  of  modifying  propellant  burn¬ 
ing  rates  Lave  achieved  limited  success.  Very  thin  metallic 
fibers  (e.g.,  aluminum,  titanium  and  zirconium)  embedded  in 
conventional  propellants  will  produce  significant  burning 
rate  increases^.  However,  these  increases  are  less  than  the 
desired  high  burning  rate.  The  overruling  objection  to  using 
metallic  fibers  is  that  no  practical  method  has  been  developed 
to  incorporate  them  into  thin  propellant  webs.  When  efforts 
are  made  to  manufacture  thin  web,  internal  burning  perforations, 
the  casting  processes  aligns  the  fibers  parallel  to  the  burn¬ 
ing  surface  so  that  they  are  not  effective  burning  rate 
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Fig.  3  Influence  of  motor  diameter  and  operating  time 
on  propellant  web  fraction. 
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enhancers.  Also  since  the  fibers  accelerate  the  burning  rate 
locally  at  a  relatively  widely  spaced  points,  uniform  burnout 
of  the  propellant  is  not  obtained. 

To  obtain  high  loading  densities,  reasonably  high  port 
area  to  throat  area  ratios,  and  high  performance  in  impulsive 
thrusters,  it  is  necessary  to  operate  at  high  chamber  pressures 
(e.g.,  greater  than  5,000  psi) .  At  these  pressures  several 
propellants  demonstrate  abrupt  changes  in  their  burning  char¬ 
acteristics,  for  example  the  burning  rate  of  some  ammonium 
perchlorate  (AP)  composite  propellants  is  known  to  rapidly 
increase  with  pressures  above  7,000  psi6/7/6.  These  rapid 
changes  in  exponent  are  associated  with  the  thermal  stress 
failure  of  oxidizer  particles  and  are  reasonably  reproducible. 
However,  the  burning  rate  exponent  greater  than  unity  (i.e., 
dlnr/dlnp  >  1)  has  discouraged  the  use  of  such  propellants. 

Another  approach  which  has  received  attention^-17  js  the 
use  of  porous  propellants .  Figure  4  shows  two  methods  for 
installing  porous  propellant  grains.  The  important  rate  con¬ 
trolling  process  in  porous  propellants  is  convective  heating 
of  the  walls  of  the  pores  rather  than  conduction  from  the 
flame  to  the  planar  burning  surface,  as  in  the  case  of  con¬ 
ventional  propellants.  Under  a  favorable  pressure  gradient, 
hot  combustion  gases  flow  intothe  pores  in  the  porous  pro¬ 
pellant  grain.  The  primary  impetus  for  studying  porous  pro¬ 
pellants  was  the  requirements  of  a  rapidly  burning  charge  at 
the  base  of  projectiles  in  traveling  charge  guns1-/16.  Several 
investigators^/16  demonstrated  high  burning  rates  using  porous 
propellants  under  the  conditions  of  dynamic  pressurization. 

But  when  these  propellants  were  used  in  traveling  charge  guns, 
they  proved  unsatisfactory  from  a  structural  standpoint11 » „ 
There  are  several  important  differences  in  the  requirements 
for  traveling  charge  guns  and  for  impulsive  v.hrusters .  The 
two  most  important  differences  are:  the  web  thickness  of  an 
impulsive  thruster  is  generally  much  less  than  the  web  thick¬ 
ness  required  for  the  traveling  charge  guns  and  the  acceler¬ 
ation  field  of  impulsive  thruster  is  generally  less  severe 
than  the  axial  and  centrifugal  accelerations  experienced  in 
the  traveling  charge  gun.  Several  methods  foi.  making  porous 
propellants  have  been  demonstrated^/16/17.  These  methods 
include  pressing  the  propellant  ingredients  to  the  desired 
density,  cementing  balls  of  propellant  together  and  then 
pressing  to  the  desired  density,  and  making  a  high  density 
felt  out  of  thin  propellant  filaments.  The  last  two  methods 
can  be  used  to  produce  propellants  with  a  wide  range  of 
porosity.  Propellants  have  been  pressed  to  within  less  than 
2%  of  their  theoretical  density.  More  details  on  the  com¬ 
bustion  and  use  of  porous  propellants  will  be  given  in  a 
report  which  is  now  being  prepared. 

One  approach  to  achieving  improved  structural  integrity 
is  to  mechanically  reinforce  the  propellant  grain.  Schultz1^ 
reported  the  results  of  a  series  of  motor  tests  in  which  thin 
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Perpendicular  alignment  of  magnetic  staples  can  be  effected 
in  a  strong  electromagnetic  field. 


a)  CASE  BONDED  GRAIN  OF  POROUS  PROPELLANT 
(REQUIRES  ELASTIC  PROPELIANT  WITH  GOOD 
PHYSICAL  PROPERTIES) 


b)  PRIMARY  PROPELLANT  IS  ROD  OF  POROUS  PROPELLANT 
AND  SECONDARY  PROPELLANT  IS  CONVENTIONAL 
PROPELLANT  BONDED  TO  CASE  WALL 


Fig.  4  Rocket  motor  configurations  that  use  porous 
propellants . 


layers  of  propellant  were  supported  by  a  thin  screen  wire  mesh. 
The  sheets  of  reinforced  propellant  were  rolled  into  scrolls . 
Tests  demonstrated  that  the  motors  could  be  rapidiy  ignited 
and  operated  at  high  pressures.  However,  structural  integrity 
of  the  reinforced  grains  continued  to  be  a  problem.  Also,  the 
wire  mesh  that  supports  the  propellant  was  largely  unburned, 
and  thus  its  presence  greatly  decreases  the  performance  of  the 
overall  propulsion  system. 

Another  approach  which  has  been  a t tempted ^ /  20  is  to  use 
propellants  with  conventional  burning  rates  and  to  obtain  the 
necessary  increase  in  burning  surface  area  by  means  of  small 
perforated  pellets  of  propellant.  The  concept  of  using  pro¬ 
pellant  pellets  to  achieve  reproducible  impulsive  thrusts  is 
well  established  in  guns.  Since  high  gas  velocities  in  the 
chamber  and  high  accelerations  tend  to  eject  the  burning  pel¬ 
lets,  in  some  exploratory  designs  the  pellets  are  retained  in 
the  motor  by  means  of  metal  cages .  The  primary  shortcoming 
of  this  approach  is  the  necessity  for  the  cage  since  the  cage 
absorbs  heat  and  its  weight  detracts  from  the  overall  perfor¬ 
mance.  The  previously  described  difficulties  of  retaining 
pellets  may  be  overcome  in  some  applications  by  taking  ad¬ 
vantage  of  high  acceleration  fields.  A  schematic  drawing  of 
such  a  rocket  motor  is  shown  in  Pig.  5.  The  propellant  is  in 
two  forms,  a  thin  layer  of  conventional  propellant  bonded  to 
the  case  wall  and  the  remainder  of  the  propellant  is  in  the 
form  of  pellets.  After  the  rocket  motor  is  ignited,  the 
charge  bonded  to  the  case  wall  burns  as  a  conventional  internal 
burning  charge  and  thus  insulates  the  motor  case.  The  pellets 
are  consumed  in  the  usual  manner.  The  acceleration  forces 
counter  the  tendency  for  the  pellets  to  be  entrained  in  the 
flow  cf  combustion  gases  and  ejected  through  the  nozzle.  One 
possible  disadvantage  of  this  approach  is  that  the  center  of 
gravity  of  the  motor  depends  upon  how  the  pellets  orient 
themselves. 

Through  out  this  discussion  we  have  alluded  to  a  less 
novel  approach  of  applying,  in  a  scaled  down  fashion,  the 
same  techniques  as  would  be  used  to  design  conventional  size 
boosters  (see  Figs.  6  and  7).  This  is  the  approach  on  which 
we  have  concentrated  c-r  efforts  during  this  study.  In  the 
sections  that  follow,  we  will  describe  this  approach.  In 
pursuing  this  approach,  it  quickly  becomes  apparent  that  con¬ 
ventional  internal  ballistics  design  techniques  are  inaccurate 
for  predicting  and  analyzing  pressure  and  thrust  versus  time 
programs  since  over  40%  of  the  operating  time  is  dominated 
by  pressurization  and  tail-off  transients.  While  we  conclude 
as  a  result  of  this  study  that  scaled  down  booster  type  designs 
are  the  proper  approach  for  the  Hit  impulsive  thruster,  we 
point  out  that  the  very  short  operating  +ime  makes  it  man¬ 
datory  for  the  designer  to  analyze  the  impulsive  thruster  by 
properly  accounting  for  the  gas  dynamics  and  burning  rate 
transients.  These  analytical  techniques  are  developed  in 
the  next  section  ai:d  then  applied  in  the  subsequent  sections. 
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Fig.  5  Example  of  impulsive  thruster  that  uses  pellets 
to  achieve  high  mass  generation  rates. 
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WEB  THICKNESS  = 0.030  IN.  {FOR  0.310  OD  CASE) 
VOLUMETRIC  LOADING  DENSITY = 0.637 


Fig<>  7  Propellant  grain  configuration  suitable  for 

low  burning  rate  propellants  (e.g.,  ~ 4  in/sec) 


DYNAMICS  OF  IMPULSIVE  THRUSTER  PERFORMANCE 


A  major  portion  of  the  study  was  devoted  to  developing 
mathematical  models  and  computer  solutions  tc  evaluate  selected 
impulsive  thruster  configurations  and  propellants.  The  dynamic 
chamber  interactions  that  must  be  accounted  for  are  indicated 
on  Figs.  1  and  8.  Figure  1  shows  the  four  major  phases  of 
motor  operation:  1)  the  ignition  phase  where  the  igniter  heats 
the  propellant  to  ignition,  2)  the  chamber  filling  phase  where 
the  gases  from  the  igniter  and  the  burning  propellant  pressur¬ 
ize  the  chamber  to  operating  conditions,  3)  the  operating  (or 
sustain)  phase  dominated  by  the  mass  generated  from  the  pro¬ 
pellant  grain,  and  4)  the  blow-down  and  tail-off  phase  after 
the  propellant  web  burns  out.  The  nozzle  is  fitted  with  a 
blow-out  disk  which  will  seal  the  nozzle  until  a  predetermined 
pressure  is  achieved. 

The  measured  performance  of  solid  rocket  motors  during 
rapid  pressure  excursions  (such  as  occur  during  rapid  ignition 
and  variable  throat  area  operation)  differs  greatly  from  pre¬ 
dictions  made  using  steady  state  burning  rate  data  and  simple 
transient  mass  balance  solutions2!-25 4  The  very  nature  of  the 
impulsive  thruster  dictates  careful  consideration  of  the  theor¬ 
etical  and  experimental  aspects  of  the  dynamic  response  of 
burning  rate,  chamber  pressure,  flame  temperature,  and  chamber 
temperature  during  both  the  rapid  pressurization  period  (as 
large  as  107  psi/sec)  and  the  high  pressure  operation  period. 
Also,  the  influence  of  these  dynamic  effects  is  very  dependent 
on  the  type  solid  propellant  being  considered.  For  example, 
the  rapid  pressure  rise  following  a  "hard"  ignition  can  pro¬ 
duce  burning  rate  overshoots  in  excess  of  fifty  percent  above 
equilibrium  conditions.  In  related  studies2^ > 25 r  we  have 
shown  that  the  seriousness  of  these  overshoots  increases  rapidly 
when  propellants  with  high  temperature  sensitivity  of  burning 
rate  are  used. 

Studies  of  dynamic  effects  have  been  limited  by  combustion 
and  gas  dynamics  models  that  are  incomplete  and  by  difficult- 
to-interpret  data.  Two  opposing  combustion  effects  cause  the 
instantaneous  burning  rate  to  differ  greatly  from  the  steady 
state  burnin„  rate  at  the  corresponding  instantaneous  pressure. 
These  are:  1)  the  nonsteady  thermal  profile  in  the  condensed 
phase  of  the  propellant  and  2)  the  out-of-phase  blowing  effect 
of  the  reactive  gases  leaving  the  burning  surface.  At  the 
lower  burning  rate  during  the  ignition  and  pressurization  phase, 
the  thermal  wave  penetrates  further  into  the  propellant  than 
at  the  operating  pressure.  As  the  burning  rate  increases,  the 
transition  to  a  thinner  thermal  wave  occurs.  Since  the  pro¬ 
pellant  is  in  effect  preheated,  the  burning  rate  is  enhanced 
while  the  overheated  surface  layer  is  being  burned  out.  The 
preheating  effect  is  only  partially  countered  by  the  excess 
blowing  effect  in  the  flame  zone,  which  tends  to  decrease  the 
heat  feedback  from  the  flame  when  the  burning  rate  exceeds 
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the  equilibrium  burning  rate.  The  interactions  resulting  from 
the  rapid  changes  in  pressure  will  also  affect  the  flame  temper¬ 
ature  and  chamber  gas  temperature,  which  will  affect  in  turn 
the  mass  discharge  rate  through  the  nozzle. 

In  previously  published  studies,  either  the  interactions 
between  the  chamber  and  propellant  combustion  were  not  considered 
or  the  burning  model  was  not  realistic  for  composite  propellants. 
In  our  development,  the  dynamic  burning  rate  law  will  be  based 
on  the  Zeldovich-Novozhilov  (Z-N)  method.  The  Z-N  method  has 
been  extensively  developed  in  the  USSR  but  has  only  recently 
been  investigated  in  the  U.S.24.  We  could  use  a  properly  for¬ 
mulated  flame  model  for  solid  propellant  combustion,  but  we 
believe  the  Z-N  method  is  most  suitable  for  *-his  application 
since  it  provides  a  more  direct  method  of  considering  complex 
propellants  in  which  the  flame  zone  parameters  are  not  known. 

The  chamber  pressure  and  mass  discharge  responses  are  calcu¬ 
lated  from  a  mathematical  model  whose  primary  components  are 
coupled  through  numerical  solutions  to  the  energy  equation  for 
the  propellant,  the  Z-N  equations,  and  the  energy  and  continuity 
equations  for  the  chamber,  which  include  the  effects  of  throat 
ablation  and  chamber  volume  charges. 


We  have  emphasized  the  dynamic  burning  and  pressurization 
of  the  impulsive  thruster  because  we  believe  that  these  are  the 
key  elements  in  evaluating  and  predicting  the  performance  of 
Hit  motor  designs  that  use  high  energy  smokeless  prooellants . * 


Analysis  of  Chamber  Flow  and  Burning  Rate  Dynamics 


The  analysis  and  assumptions  are  discussed  as  they  apply 
to  impulsive  thrusters  with  internal  burning,  perforated  grains, 
e.g.,  pressures  up  to  20,000  psi,  burning  rates  in  the  range 
of  2  to  20  in/sec  (over  the  range  of  pressures) ,  and  dp/dt  up 
to  5  x  10 7  psi/sec.  Because  of  the  wide  range  of  events  and 
special  purpose  propellants  that  may  be  considered,  the  assump¬ 
tions  used  in  the  analysis  should  be  reviewed  for  each  situation- 
considered.  As  indicated  in  Fig.  8,  we  are  considering  two 
interacting  sets  of  dynamic  processes:  1)  the  flew  of  gases  in 
the  chamber  and  2)  the  propellant  combustion.  The  assumptions 
and  equation  development  for  these  processes  are  treated 
separately. 


The  assumptions  for  the  propellant  are: 


1)  The  rate  processes  in  the  gas  phase  and  in  the  surface 
reaction  zone  can  be  considered  quasi-steady  in  the 
sense  that  their  characteristic  times  are  short  compared 
to  the  transition  time  of  the  pressure  change.  No  such 
limitation  is  necessary  for  the  thermal  wave  in  the  con¬ 
densed  phase. 

2)  No  kinetic  heat  release  occurs  in  the  condensed  phase 
below  the  surface  reaction  zone.  Although  the  surface 

*  As  will  be  pointed  out  in  subsequent  sections,  a  properly  de¬ 
signed  motor  will  operate  in  a  domair.  where  these  affects  are 
not  detrimental.  Thus,  one  of  the  objectives  is  to  areatlv 
reduce  the  very  effects  we  are  emphasizing. 


reactions  occur  in  a  zone  of  finite  thickness,  the  zone 
is  sufficiently  thin  that  it  can  be  considered  as  quasi¬ 
steady. 

3)  The  condensed  phase  of  the  propellant  is  accurately  repre¬ 
sented  as  being  homogeneous  and  isotropic. 

4)  Propellant  combustion  zones  are  not  influenced  by  axial 
position  and  external  forces  such  as  shear  forces  from 
the  flowing  chamber  gases  and  acceleration. 

The  assumptions  for  the  chamber  are: 

5)  Effects  of  axial  variations  of  chamber  temperature  and 
pressure  can  be  accounted  for  by  defining  a  mean  burning 
rate  point  between  the  head  end  and  the  nozzle  end  of 
the  chamber. 

6}  The  chamber  gases  can  be  treated  as  perfect  gases. 

?)  Heat  losses  to  the  inert  parts  are  negligible. 

8)  The  specific  heats  and  average  molecular  weights  for  the 
chamber  gases,  flame  and  igniter  gases  are  equal. 

For  the  conditions  considered  in  this  study,  the  justifi¬ 
cations  for  assumption  1  closely  follow  the  arguments  pre¬ 
sented  in  Ref.  24.  Thus  following  Ref.  24,  the  characteristic 
times  for  the  condensed  phase,  surface  zone,  and  gas  phase  are 

t c  =  uc/r2  55  10“4  sec  for  r  =  2.0  in/sec  (1) 

~  10" 6  sec  for  r  =  20.0  in/sec 


xs  =  R  Ts/E  xc  ~  0.1  Tc  <  i0~5  sec  (2) 

Tg  =  [Ag  ccpg/(Ar  cgpc)3  t c  -  0.01  t c  <  10'b  sec  (3) 


Since  the  times  of  the  monotonic  portion  of  the  pressure  excur 
sion  considered  in  this  paper  are  approximately  0.001  sec., 
the  conditions  of  assumption  1  are  satisfied. 

The  arguments  for  assumptions  2,  3  and  4  are  the  same  as 
presented  in  Ref.  24. 

The  condition  of  negligible  axial  temperature  variations 


cf  assumption  5  is  more  nearly  justified  for  an  internal¬ 
burning  propellant  grain  with  a  low  L*  than  for  an  end-burning 
propellant  grain  wit>  a  large  L*.  In  the  former  case  which 
is  the  situation  in  this  study,  the  combustion  gases  leaving 
the  burning  propellant  surface  are  distributed  along  the  length 
of  the  rocket  motor  md  have  radial  velocity  components  (normal 
to  the  axis)  which  promote  rapid  mixingf  of  these  gases  with 
the  gases  flowing  from  the  head  end  of  the  motor.  Also,  in 
low  L*  combustors  the  stay  times  of  the  combustion  gases  in  the 
chamber  are  considerably  less  than  the  transition  time  for 
the  pressure  change,  and  as  a  result  large  axial  thermal  gra¬ 
dients  can  not  develop.  For  example,  the  characteristic  stay 
time  of  the  datum  case  chamber. 


ch 


=  V„/A+ 


(4) 


is  0.0003  sec  while  the  monotonic  portion  of  the  pressure  transi¬ 
tion  is  approximately  0.001  sec.  Also,  axial  temperature  vari¬ 
ations  are  reduced  by  the  mixing  produced  by  turbulence , 
free  convection  currents,  and  vortices.  It  should  be  noted 
that  the  stipulation  that  L*  be  small  is  not  a  restriction  to 
the  analysis  since  dynamic  effects  are  prominent  motors  with 
relatively  low  L*  values.  The  limitation  imposed  by  assump¬ 
tion  5  is  that  the  propellant  can  not  be  susceptible  to  erosive 
burning.  Because  we  are  considering  very  high  burning  rate 
propellants ,  the  analyses  of  Willoughby20  and  Saderholm2'  indicate 
that  erosive  burning  will  not  occur.  Saderholm’ s  criterion, 
which  has  been  successfully  applied  to  a  number  of  high  burning 
rate  composite  propellants,  indicates  that  for  lach  numbers  less 
than  0.3  and  chamber  pressures  above  10,000  ps  erosive  burning 
will  not  occur  if  the  basic  burning  rate  is  in  excess  of  3.C 
in/sec. 

Assumption  6  is  widely  employed  in  the  combustion  litera¬ 
ture  for  the  conditions  considered  in  this  study. 

Assumption  7  is  consistent  with  the  small  amount  of  exposed 
motor  hardware  in  properly  designed  high  performance,  internal 
burning  systems.  Also,  the  small  amounts  of  refractory  insula¬ 
tions  that  are  used  absorb  very  little  heat.  This  is  demon¬ 
strated  by  the  delivered  specific  impulses  of  well  designed, 
high  operating  pressure  motors  having  delivered  specific  im¬ 
pulses  that  are  within  5%  of  the  theoretical  specific  impulses. 


t  These  flow  patterns  have  recently  been  photographed  in  a 
high  L/D,  low  L*  window  burner  which  is  being  developed  at 
Princeton  to  study  ignition  transients.  The  rapid  mixing  of 
the  gas  is  evident. 


The  conditions  of  Assumption  8  were  tested  for  the  datum 
case  usina  the  calculated  equilibrium  thermochemical  results 
of  Ref.  28.  The  values  of  specific  heat  and  molecular  weight 
varied  less  than  5%  over  the  500  F  range  of  temperature. 

Functional  Relationships  for  the  Chamber 

The  previously  discussed  pressurization  sequence  is  shown 
in  Fig.  1  and  chamber  interactions  are  indicated  in  Fig.  8. 

The  initial  conditions  are 


r  =  0 


?ch  pinit 

T  =  T . 
ch  imc 


(5) 


At  =  0 


i.e. , 


nozzle  closure 


in  place 


The  nozzle  is  blocked  by  a  nozzle  closure  which  completely 
seals  the  throat  until  a  predetermined  chamber  pressure,  Pbiow< 
is  attained.  Then  the  throat  area  is  allowed  to  increase 
linearly  in  a  prescribed  time  interval  for  the  nozzle  closure 
to  clear  the  throat,  Atayea*  During  motor  operation  the  throat 
diameter  is  permitted  to  increase  due  to  erosion. 


AD.  =  2 


dt 


(6) 


then  in  turn 


At  “  lDt,init  +  “t> 


(7i 


The  mass  continuity  equation  for  the  free  volume  in  the 
motor  chamber  (see  chamber  control  volume  in  Fig.  8)  is 


d<Pchvch>/dt  +  =  mb  +  mig 


(8) 


where  for  a  perfect  gas 


”n  =  PchAt  g/c* 


(9) 
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=  pcAb(t)  r(t) 


(10) 


During  the  ignition  phase,  the  burning  surface  area  increases 
as  the  flame  spreads  along  the  propellant  surface;  following 
full  ignition  the  burning  surface  area  is  a  prescribed  func¬ 
tion  of  the  grain  geometry  in  terms  of  distance  burning.  Thus, 


for  J  vfs  dt  /  L 


prop 


<  1 


(11) 


(12) 


By  expressing  PChvch  terms  °f  the  perfect  gas  law, 
the  following  relation  for  dpcb/dt  is  obtained  from  Eq.  8 


dpch/dt  =  -(Pch/vch>dvch/dt 

+  (YR/Vch)  (--mnTch  +  mbTf  +  migTig) 

The  energy  equation  for  the  free  volume  in  the  motor 
chamber  is 

d(pchvchTch)/dt  +  mnTch  ■  «bTf  ‘  ™igTig 


(l/cp)d(pchVch)/dt  =  0  . 


(14) 


_  Ol  _ 


Multiplying  the  mass  continuity  equation,  Eq.  8,  by  Tch  and 
subtracting  it  from  the  energy  equation,  Eq.  14,  multiplied 
by  y  yields 


dWdt  =  -(lRn/pchVch)Tch(Y“i) 


+  (mn/pchVch)  (YTf  -  Tch  +  YTigmig/mb) 


(15) 


The  chamber  volume  increases  as  the  propellant  is  burned: 


dVch/dt  =  rAb 


(16) 


The  several  regions  in  the  motor  which  were  coupled 
through  one  dimensional  compressible  flow  solutions  are  shown 
in  Fig.  9.  The  solution  for  pc*  and  T  ,  applies  to  the 
mean  flow  condition  for  the  chamber  and  satisfies  A  -  VC/  an£ 
Aoort  mean*  From  these  mean  conditions,  the  stagnation  con- 
dltiofis''at  the  head  end  pb,  Aport  ,h'  ' Tb,h  and  rb  and  the 
nozzle  end  conditions  p.. -'p0  ,  Mn'and  A'*  n,  tb  nand  r 
are  calculatLoned  using  the  1-D  dimension  compressions  flow11 
(with  mass  addition)  equations.  The  thrust  calculations 
are  based  on  the  stagnation  pressure  at  the  nozzle,  p^n  . 

Functional  Relationships  for  Propellant  Combustion  Transients 

In  this  development,  the  dynamic  burning  rate  law  is 
based  on  the  Zeldovich-Novozhilov  (Z-N)  method.  The  Z-N 
method  offers  important  advantages  when  considering  the  burn¬ 
ing  rate  transients  of  propellants  whose  burning  rate  mec¬ 
hanisms  are  not  understood  (e.g.,  highly  catalyzed  propellants, 
propellants  with  complex  r  vs  p  relationships,  and  propel¬ 
lants  with  high  percentages  of  metal  fuels) .  It  is  important 
to  note  that  the  valid’ ty  of  applying  the  Z-N  method  to  a 
particular  combustion  situation  can  only  be  established  by 
flame  zone  studies  that  determine  the  applicability  of  the 
assumptions.  Flame  zone  studies  revealed  that  the  Z-N  method 
is  valid  for  the  range  of  conditions  being  considered  for  Hit. 

We  emphasized  in  Ref.  24  that  the  usefulness  of  the  Z-N 
method  was  limited  because  of  the  lack  of  Ts(P/'I’n)  data.  As 
we  planned  experiments  to  develop  this  data,  it  became  appar¬ 
ent  that  the  concept  of  a  surface  temperature  in  the  sense 
that  Zeidovich  meant  it  is  misleading.  The  steady  -state  data 
which  is  used  in  this  development  of  the  dynamic  burning  rela¬ 
tionship  is  correlated  by  the  pyrolysis  law 
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r  =  As  exp(-Ss/RTif )  (16) 

which  is  based  on  the  temperature  at  the  interface  between 
the  very  thin,  quasi-steacy  surface  reaction  zone  (see 
assumption  2  and  Fig,  8)  and  the  nonreacting  condensed  phase. 
(It  is  important  to  note  that  under  dynamic  conditions  T ■ 
is  not  a  single  valued  function  of  pressure  but  is  depemMnt 
on  the  pressure  and  burning  rate  history.  This  is  a  departure 
from  previous  applications  of  the  Zeldovich  method  (and  later 
the  Z-N  method)  in  which  a  surface  temperature,  Ts,  (i.e., 
the  temperature  at  the  gas/burning  surface  interface)  was 
used  in  the  pyrolysis  law.  Once  it  is  realized  that  the 
surface  thermal  zone  is  not  a  well  defined  zone,  it  is  more 
meaningful  (and  less  arbitrary)  to  base  the  pyrolysis  law 
on  rather  than  T  . 

The  energy  equation  for  the  condensed  phase  (-»  <  xc  <  0) 
shown  in  Fig.  8  has  the  following  eigenvalue  dependence  on 

T  ( 0 ,  t )  =  T  (t) 
c  if 


pcccOTc/3t  +  r(Tif)  8Tc/3x]  =  Ac  32Tc/3x2  (17) 


where  under  conditions  of  full  ignition  the  burning  rate  r 
is  related  to  the  interface  temperature  through  Eq.  16. 

The  initial  condition  for  Eq.  17  is. 


T(x)  =  TQ  at  t  =  0  (18) 

The  first  boundary  condition  is 


3T/3x  =0  as  x  -*•  -»  (19) 


The  second  boundary  condition  is  a  series  of  conditions: 

1)  heatup  to  gasification 

Ac(3T/3r)c,if  =  qign  for  0+  <  t  <  ty  (20) 

2)  gasification  prior  to  establishment  of  flame 

=  Sign  -  rpAhs  for  fcv  <  ±  fcz  eld 


(21) 


-  9/1 


ENERGY  BALANCE  AT  SURFACE: 


rpccc(Tif  "V  =  'c«  =rccQs  -"‘(‘f.if 


STEADY  STATE  EXPERIMENTAL  RESULTS: 


TRANSFORMATION  TO  r  -  0  COORDINATES : 


0  (k/pc) _  {Tif~T0) 
c 


^MEASURED 
p_T  POINTS 

UL* 


Fig.  10  Schematic  drawing  showing  relationship  between 
experimental  data  and  the  required  heat  feed¬ 
back  function  for  the  dynamic  burning  rate  model. 
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3}  combined  heating  from  ignition  stimulus  and  from  gas 
phase  and  surface  reactions 

=  g.  +  X  <P  (r,p(t) )  for 
ign  c 


t '  <  t  <  t.  (22) 

2eld  -  ign,otf 


4)  burning  during  conditions  of  full  ignition 


=  A  P (r ,p (t) )  for  t  >  t.  , 
c  ign'ofr 


(23) 

where  <J>  is  the  2-N  heat  feedback  function  developed  in  Ref. 

24  except  that  it  is  based  on  rather  than  T  . 

Figure  10  contains  a  concise  outline  of  how  the  heat 
feedback  function  i  is  developed. 

The  following  relationships  are  available  in  principle 
from  laboratory  experiments 

r°  =  r°(T0,p)  (24) 

and 

T®  =  T®  (p»r)  (25) 

if  if 

The  information  of  Eq.  24  can  be  obtained  by  routine  burning 
rate  experiments  (for  a  limited  temperature  range) .  Experi¬ 
mental  measurements  of  the  dependence  of  on  p  and  r 
are  much  more  difficult  and  have  yet  to  be*accomplished  with 
good  accuracy  over  the  desired  ranges  of  r  and  p.  The 
Z-N  method  uses  steady-state  burning  rate  data  in  conjunction 
with  the  following  expression  for  the  steady-state  energy 
balance  in  the  condensed  phase 


,cd  =  r°cccc(T?f 


V 


(26) 


As  indicated  in  Fig.  10,  rfT-*:)  and  r(p,T  )  data  can  be  combined 
vith  Eq.  26  to  obtain  the  following  functional  form,  which  re¬ 
places  the  corresponding  function  which  is  normally  derived 
from  flame  model. 


_  A  O  —O  \ 

-  vU  /P  ; 


P[r(t)  ,p(t)J  . 


(27) 


Tnis  relation,  which  is  obtainable  from  steady-state  experi¬ 
ments,  is  the  key  element  of  the  2-N  method,  since  it  expresses 
the  heat  feedback  to  the  condensed  phase  in  terms  of  two  para¬ 
meters,  pressuz'e  (p)  which  is  externally  imposed,  and  burning 
rate  (r)  which  is  the  important  explicit  result  desired  from 
the  solution.  Assumption  1  is  the  justification  for  applying 
the  steady-state  result  of  Eq.  28  to  transient  conditions. 


At  a  given  pressure,  there  can  be  a  wide  range  of  burning 
rates,  r:  1)  during  steady-state  conditions,  r  is  varied 


Based  on  the  arguments  xn  Ref.  Z4,  it  is  maintained  tnat  for  a 
given  r  and  p  (regardless  of  how  the  value  of  r  is  ob¬ 
tained)  the  value  of  $(r,p)  is  the  same  for  both  the  steady- 
state  and  transient  conditions.  It  is  important  to  realize 
that  even  though  Eq.  27  is  obtainable  from  steady-slate  data, 
transient  burning  rate  calculations  can  be  performed  only  if 
the  transients  in  the  condensed  phase  are  calculated,  (i.e., 
•£(t))  by  solving  the  transient  heat  condition  equation  (Ec.  16) 
with  its  boundary  conditions.  Thus,  the  chamber  pressure,  p(t), 
is  one  of  the  dependent  variables  to  be  obtained  from  the 
simultaneous  solution  of  the  energy  and  continuity  equations 
for  the  chamber. 


An  expression  for  the  dynamic  flame  temperature  may  be 
obtained  under  the  following  conditions:  1)  the  adiabatic 
product  composition  in  the  flame  zone  and  the  heats  of  re¬ 
action  in  the  surface  reaction  zone  and  in  the  flame  zone  are 
not  affected  by  the  dynamic  conditions,  and  2)  the  specific 
heats  in  the  surface  reaction  zone  and  the  flame  are  equal. 

It  follows  from  b  =  (r/ac)  (T^c-Tq)  that  the  energy  balance 
for  both  steady  and  unsteady  stages  is: 


roc£hf  =  XcC  +  roccf(Tf-Tif). 


Therefore  by  considering  a  reference  condition. 


Dccf  M, rJ  ref 


-  | T.  „ 

4  1  *■  . 


S  fc  I 

-  I  *  C  ; 


The  result  of  the  foregoing  analysis  of  the  chamber  and 
propellant  burning  conditions  is  a  coupled  set  of  three  non¬ 
linear,  first  order  ordinary  differential  equations  (Eq.  13, 
15  and  16)  for  Pqq/  Tck,  and  vch  an°  one  nonlinear,  second 
order  partial  dir ierer. rial  equation,  {Eq.  17)  for  T(x). 


Equation  17  with  its  boundary  conditions  was  solved  using  the 
methods  of  explicit  finite  differences.  Equations  13,  15  and 
16  were  integrated  by  the  fourth  order  Runge-Kutta  technique 
using  a  variable  time  step.  For  special  limiting  cases,  the 
numerical  results  agreed  with  closed  form  solutions  for  cham¬ 
ber  venting,  adiabatic  compression,  etc. 


Example  of  Measured  Dynamic  Burning  Rates 


In  a  related  study,  transient  pressure  data  (from  which 
the  dynamic  burning  rates  can  be  determined)  were  obtained 
by  suddenly  reducing  the  nozzle  throat  area  of  a  motor  with 
a  rcd-and-tube  propellant  charge  designed  to  have  a  very  small 
free  volume  in  order  to  have  a  rapid  pressure  rise2-5.  (This 
is  referred  to  as  the  Low  L*  Combustor.)  After  ignition,  the 
motor  is  allowed  to  reach  equilibrium  at  low  pressure,  and 
then  the  throat  area  is  rapidly  reduced  by  insertion  of  a 
pintle  driven  by  a  pneumatic  cylinder.  The  chamber  pressure 
rises  sharply  to  a  new  equilibrium  level.  This  method  of  ob¬ 
taining  rapid  pressure  changes  and  dynamic  burning  rate  data 
gives  a  more  direct  result  than  previously  attempted  methods. 
For  example,  Wooldridge  and  Marxman22  used  propellant  powders 
(as  a  secondary  charge)  to  pulse  their  combustor,  as  well  as 
a  rapid  decrease  of  the  threat  area.  Establishing  the  dynamic 
burning  rate  from  such  an  experiment  is  complicated  by  the 
impulsive  changes  brought  on  by  the  secondary  charge.  Similar¬ 
ly,  obtaining  dynamic  burning  rate  information  from  T-burner 
experiments  is  complicated  by  approximations  cf  the  losses  in 
the  T-burner25.  Also,  it  is  not  feasible  to  deduce  dynamic 
burning  rate  information  from  the  rapid  pressure  rise  follow¬ 
ing  ignition  since  the  complexities  of  ignition  processes 
preclude  isolating  the  desired  dynamic  burning  effects. 


Figure  11  shows  a  rapid  pressurization  and  overshoot  (in 
the  case  of  test  2)  resulting  from  the  dynamic  super-rate  ef¬ 
fects  for  a  sudden  throat  area  decrease  (in  0.001  sec).  Al¬ 
though  test  1  also  exhibits  a  strong  super-rate  effect,  the 
overshoot  of  pressure  effect  is  muffled  by  the  somewhat  larger 
L*  associated  with  the  smaller  final  A*..  The  propellant  was 
70%  AP,  29%  PBAA,  1%  CuO  &  (^203.  Other  motors  have  been 
tested  over  a  range  cf  chamber  conditions.  The  instantaneous 
burning  rates  and  chamber  temperatures  were  deduced  from  the 
measured  p  vs  t  by  means  of  a  solution  to  the  energy  and  con¬ 
tinuity  equations  for  the  chamber.  This  solution  is  for  a 
much  simpiier  situation  than  the  one  previously  described  in 
this  section  since  flame  temperature  (not  the  chamber  temper¬ 
ature)  is  assumed  to  be  constant  and  the  partial  differential 
equation  governing  heat  flow  to  the  condensed  phase  is  not 
solved.  Figure  12  shows  the  dynamic  burning  rate,  r/r„a, 
determined  from  tesc  1.  In  Fig.  13,  r/r  is  displayed  "in 
terms  cf,fche  dimensionless  p  parameter  which  was  used  bv 
Von  Elbe  .  The  dynamic  burning  rates  of  Fig.  13  demonstrate 
a  hysteresis-like  character  rather  than  a  simple  straight- 
line  relationship  between  the  dimensionless  p  and  dynamic 


Comparison  of  calculated  and  experimental  pressure-time  trace 
for  two  A.  values  (example  dynamic  burning  and  pressurization 
effects). 


TIME,  t,  MSEC 


Fig.  13  Experimentally  determined  burning  rate 
correlated  with  a  dimensionless  p. 
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burning  rate,  as  given  by  some  of  the  previously  developed 
dynamic  burning  rate  equations. 

Rocket  Motor  Performance  Parameters 

The  basic  inputs  to  the  analysis  are  desired  average 
thrust,  total  impulse,  propellant  properties,  and  nozzle 
characteristics.  These  inputs  were  selected  since  they  are 
frequently  specified  as  requirements  to  achieve  a  particular 
mission.  The  motor  performance  equations  are  presented  to 
avoid  confusion  concerning  how  the  motor  parameters  were 
applied. 

The  steady-state  burning  rate  along  the  grain  is  ex¬ 
pressed  as  a  function  of  static  pressure,  p,  and  grain 
temperature ,  Tq / 

r  =  a  pn  exp[op(Trof-T0) ]  (30) 

The  values  a,  n  and  o  are  measured  properties  of  the  pro¬ 
pellant.  The  factor 


exp[ap(Tref-T0)J 

accounts  for  the  changes  in  steady-state  burning  rate  caused 
by  the  departure  of  ;he  grain  temperature,  Tq,  from  the  ref¬ 
erence  temperature  Trej. 

A  compatible  set  oc  combustion  gas  properties  is  estab¬ 
lished  by  starting  with  calculated  combustion  gas  properties 
(obtained  from  thermochemistry  computer  programs  such  as 
Ref.  28)  and  then  adjusting  the  calculated  properties  to  ob¬ 
tain  agreement  with  measured  mass  flow  rates  and  thrusts. 

The  measured  value  of  c*  that  is  normally  used  to  calculate 
the  mass  discharge  is  established  from  motor  firing  data  by 
the  following  equation 


9  At  1  pdt 
cmeas  =  : 

wprop 


(31) 


The  quantity  c*  is  used  to  relate  instantaneous  mass  flow  rate, 

m  to  throat  area  A  and  chamber  pr.  ssure  p 
n  t  on 


(32) 
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From  isentropic  one-dimensional  compressible  flow  theory, 
mass  discharge  rate  can  be  expressed  in  terms  of  the  gas 
properties  as  follows. 
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Ideally,  c*  is  defined  in  terms  of  the  chamber  gas  properties 
as 
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To  have  a  compatible  set  of  relationships,  the  interrelation 
between  c* 


and 


_  _  -  T  .  4.  that  is  established  by  Eq<  34  is 

maintained  throughout0 “analysis .  The  measured  value  of 
is  input  and  T^  act  is  calculated  from  Eq.35 
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It  is  interesting  to  note  that  the  measured  c*  values  corrected 
to  ideal  standard  conditions  are  within  94  to  97  percent  of  the 
theoretical  c*  calculated  from  thermochemical  programs  such  as 
Ref.  2s* 

The  ratio  of  specific  heats,  y,  is  an  average  effective 
value  for  the  motor  and  nozzle.  When  possible  this  '  alue  was 
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and 


X  ={l+cos  a  )/2 
n 


(40) 


The  value  of  motor  coefficient,  C^,  used  in  the  equation  is 
for  nozzle-end  stagnation  pressure  and  should  not  be  confused 
with  values  that  were  determined  from  the  relationship  between 
head-end  pressure  and  measured  thrust. 

Total  Impulse  of  Motor  Rotating  in  Plane  of  Desired 
Thrust  Vector 

With  the  specification  of  a  thrust  time  trace  directed 
radially  from  the  axis  of  rotation  of  the  vehicle,  it  is 
desirable  to  specify  the  time  at  which  ignition  should  occur, 
such  that  maximum  thrust  is  delivered  in  a  given  direction 
(see  Fig.  14).  In  this  discussion  the  vehicle  rotates  with 
constant  angular  velocity  as,  clockwise,  and  thrust  is  to  be 
maximized  in  the  +y  direction.  The  given  thrust  time  trace 
has  total  impulse  (when  stationary)  of 
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fin  Fdt 


(41) 


blow 

The  radial  position  of  the  vehicle  is  described  by: 

0=a)t+e^ 


(42) 


where  6  is  the  angle  formed  by  the  thrust  vector  and  the  +y 
axis.  It  is  assumed  that  as  is  such  that  motor  burn-out  takas 
place  in  a  fraction  of  one  revolution. 

A  rotational  total  impulse  is  defined  Dy: 


I  =(F  ,F  )  = 
r  x  y 
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and  F  is  to  be  maximized,  while  F  should  be  small. 

y  ^ 

A 

A  Substituting  for  6  in  I  ,  and  solving  for  the  magnitutde 
of  Ir, 
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Since  j  I  j  is  constant  for  a  given  as  and  always  less 
than  Im,  it  may  be  regarded  as  the  maximum  effective  thrust 
attainable  with  a  giver,  oj. 


AND  9  =  C  COINCIDE 
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Fig.  14  Schematic  drawing  defining  nomenclature  used  in 

analysis  of  ignition  delay  errors  on  total  impulse 
vector. 
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Define  the  centroid  time  of  the  thrust  versus  time 

program  ,  t  ,  such  that: 
cen 

f" cen  j*  fin 

J  Fdt=J  Fdt  (45) 

V  blow  tcen 


A  series  of  calculations  were  performed  to  demon¬ 
strate  the  influence  ignition  delay  errors  on  the  magni¬ 
tude  and  direction  of  the  thrust  vector.  The  ignition 
delay  error  is  defined  as  the  time  interval  between  the 
instant  the  nozzle  centerline  is  aligned  wich  the  de¬ 
sired  direction  of  the  total  impulse  and  the  centroid 
time,  t  ,  defined  by  Eq.  45.  Figure  15  is  the  thrust 
versus  ??me  program  used  in  the  calculations.  Figure  16 
shows  the  influence  of  the  ignition  delay  error  on  the 
thrust  vector.  Note  that  an  ignition  delay  error  of 
0.00125  sec.  (the  3a  value  in  Table  1)  results  in  the 
thrust  vector  being  misaligned  by  11  degrees.  The  thrust 
smear  factor,  defined  as 


F  =  I  /I 
s  r  T 


has  a  maximum  value  of  0.95  when  ignition  delay  is  zero 
and  falls  off  slightly  to  0.93  when  the  ignition  error  is 
0.00125  sec. 


(46) 


In  the  section  entitled  "Performance  Analyzed  in  Terms 
of  F  vs^t",  the  effects  of  motor  and  propellant  variables  on 
F  and  Tj.  will  be  considered  for  a  particular  design. 


*3  O 

v/  u 


HIGH  BURNING  RATE  PROPELLANTS 


General  Consid6i at  t  _ s 

A  wide  range  of  high  burning  c^te  propel la- - *  -ere 
investigated.  The  most  explicit  cor.-. i;  amt  ;  •-  -  pro¬ 
pellant  is  that  it  have  a  burning  rare  of  at  least  6  in/sec 
at  10,000  psia.  As  indicated  if-  T-ahlj  2  se  eral  propellant 
systems  can  meet  this  requirement-  The  other  items  which 
influence  the  propellant  selection  azv  'ore  subtle  and 
must  be  evaluated  ;n  terms  of  system  requirements  and  motor 
trade-offs . 

As  with  most  volume  limited  systems,  we  are  seeking 
a  propellant  with  a  maximum  IgDo  product.  This  maximum 
must  be  obtained  without  using*  elemental  rnetal  fuel  addi¬ 
tives  such  as  aluminum.  Independent  of  exhaust  plume 
considerations,  elemental  metal  additives  are  excluded  be¬ 
cause  they  form  metal/metal-oxide  agglomerates  that  are 
likely  to  restrict  the  throat  and  to  present  serious  two 
phase  flow  problems  in  the  very  small,  center  vented  cham¬ 
ber.  Generally  the  dual  requirements  of  high  Isp  and 
burning  rate  are  compatible  since  high  oxidizer  loadings 
increase  both  burning  rate  and  specific  impulse.  However, 
the  upper  limit  of  propellants  1  and  2  in  Table  2  may  de¬ 
pend  on  processing  and  casting  limitations  resulting  from 
the  high  specific  surface  of  fine  ammonium  perchlorate  (AP) • 
Greater  flexibility  in  burning  rate  and  more  nearly  opti¬ 
mum  ammonium  perchlorate  packing  fractions  can  be  obtained 
by  using  an  iron  compound  burning  rate  catalysts  {e.g., 
propellant  No.  2). 

A  primary  consideration  in  the  propellant  selection 
will  be  the  reproducibility  and  predictability  of  the  cent¬ 
roid  of  the  thrust  versus  time  trace.  Reproducibility  prob¬ 
lems  can  originate  during  the  flame  spreading  phase.  Based 
on  our  present  understanding,  a  propellant  with  a  very  rapid 
flame  spreading  rate  is  necessary  since  the  rapid  flame 
spreading  characteristic  will  dominate  the  ignition  process 
and  tend  to  overcome  effects  of  nonuniform  heating  by  the 
igniter.  Also  if  the  flame  spread  period  is  small  (i.e., 
less  than  0.0002  sec.)  relatively  large  percentage  deviations 
in  the  flame  spread  rate  will  not  produce  large  deviations  in 
the  ignition  time  interval.  Another  advantage  of  an  easily 
ignited  propellant  is  that  the  igniter  weight  will  be  small. 
Indeed,  it  is  expected  that  by  proper  propellant  selection 
the  Hit  impulsive  thruster  can  be  ignited  by  a  squib  with  a 
charge  on  the  order  of  0.1  gm.  Propellant  No.  2  is  known  to 
have  rapid  flame  spread  characteristic  and  to  require  a  rela¬ 
tively  low  ignition  energy. 

The  exha’-st  plume  considerations  on  propellant  selection 
involve  several  overlaping  objectives  and  will  be  discussed 
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in  a  separate  section  of  the  report. 

Survey  of  High  Burning  Rate  Propellants 

A  wide  range  of  solid  propellants  have  been  examined  in 
terms  of  the  Hit  requirements.  On  17  May  1571,  M.  Summerf ield 
and  L.  H.  Caveny  of  Princeton  visited  the  Feltman  Laboratories 
of  Picatinny  Arsenal  and  discussed  the  propellant  selection 
with  J.  P.  Picard,  C.  Lenchit^,  R.  Wetton  and  R.  P.  Baumann. 

On  11  May  1971,  L.  H.  Caveny  visited  the  Naval  Ordnance  Station 
at  Indian  Head,  Maryland  and  discussed  the  Hit  propellant 
selection  and  internal  ballistics  with  A.  Camp  and  J.  H. 

Wiegand.  Since  both  the  Feltman  Laboratories  and  the  Naval 
Ordnance  Station  are  primarily  concerned  with  double-base 
propellants,  the  discussions  centered  on  what  can  be  done  to 
increase  the  burning  rate  and  performance  of  double -base 
propellants.  It  was  generally  agreed  that  the  standard  M-7 
propellant  represents  an  upper  limit  of  burning  rate  combined 
with  good  specific  impulse.  As  a  result  of  these  discussions 
a  series  of  M-7  strands  was  burned  at  the  Feltman  Laboratories. 
At  10,000  psi,  the  strand  burning  is  4.0  in/sec.  This  is 
probably  the  upper  limit  of  burning  rate  since  M-7  and  M-9 
burning  rates  in  motors  are  lower  than  strand  burning  rates. 

Thus  the  Kit  requirements  can  be  achieved  only  by  using  very 
thin  web  propellant  configurations ,  such  as  a  double  web 
internal  burning,  case-bonded  propellant  grain  with  a  single 
web  thickness  of  0.030  inches. 

On  31  March  1971,  L.  K.  Caver.y  visited  the  Huntsville 
Division  of  the  Thiokol  Chemical  Corporation  and  discussed 
high  burning  rate  AP  composite  propellants  with  W.  E.  Hunter. 

D.  A.  Flanagan  and  J.  0.  Hightower.  Much  of  the  Thiokol  inform¬ 
ation  is  included  in  Ref.51  In  response  to  letters  requesting 
specific  propellant  information  the  Atlantic  Research 
Corporation3*  and  Alleghany  Ballistics  Laboratory  of  Hercules 
Incorporated3  3  provided  information  cn  several  propellant 
types , 

Because  of  the  many  considerations,  in  our  opinion  it  was 
not  advisable  to  eliminate  at  this  stage  of  the  evaluation 
several  of  the  propellants  that  have  questionable  character¬ 
istics.  The  five  propellants  (which  we  believed  to  be 
candidates)  are  listed  in  Table  3  which  also  contains  a  list 
of  properties  which  are  needed  to  make  a  preliminary  evalua¬ 
tion  of  each  propellant.  All  of  the  properties  in  Table  3 
are  input  properties  to  our  Transient  Gas  Dynamics  and 
Burning  Rate  Model  and  Motor  Components  Model. 

Propellant  Types 

The  propellant  to  meet  the  requirements  of  Table  i  should 
have  a  burning  rate  of  at  least  six  inches  par  second  at 
10,000  psi  and  have  a  high  specific 
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impulse  without  the  use  of  metal  fuels,  since  there  is  no 
ideal  system  in  terms  of  availability/  high  specific  impulse, 
high  density,  and  low  temperature  sensitivity,  it  is  necessary 
to  consider  various  trade-offs  among  these  parameters  to 
determine  the  most  suitable  propellant.  A  bibliography  of 
modern  high  burning  rate  propellants  is  given  in  Appendix  1. 

Some  propellants  based  on  boron  compounds  have  high 
burning  rates,  high  specific  impulses,  and  excellent  tempera¬ 
ture  sensitivity  characteristics.  Tests  of  these  propel¬ 

lants  have  demonstrated  that  the  castmq  properties,  the 
bonding  capability,  and  the  physical  strength  are  adequate 
for  use  in  the  Hit  system.  It  should  be  noted  that  these 
propellants  will  have  a  high  percentage  of  boron  oxides  in 
their  exhaust.  A  potential  advantage  of  the  boron  compound 
propellants  is  that  the  boron  oxides  in  the  chamber  will 
assist  in  dissipating  some  types  of  combustion  instability. 

It  is  expected  that  these  oxides  are  in  the  sub-micrcn  size 
range.  However,  these  oxides  will  increase  the  cool  down 
time  of  the  exhaust  plume- 

NF  propellants  have  been  actively  considered  in  studies 
to  obtain  a  low  smoke  propellant.  These  propellants  contain 
no  meta]  oxides  in  the  exhaust  and  have  high  specific  impulses 
over  a  satisfactory  range  of  burning  rates.  A  difficulty  with 
the  NF  propellants  is  the  high  temperature  sensitivity  of 
the  burning  rate  ,  as  great  as  0.40%/K.  Based  cn  our  recent 
studies  of  pressure  overshoots  resulting  from  dynamic  burning 
rate21*'25  we  expect  that  propellants  with  temperature  sensitiv¬ 
ity  values  in  this  range  will  have  serious  ignition  problems. 

One  of  the  problems  experienced  in  obtaining  a  truly  smokeless 
NF  propellant  was  the  condensation  of  exhaust  products  at  sea 
level  conditions.  Under  the  vacuum  conditions  cf  Hit,  exhaust 
condensation  is  unlikely.  Thus,  the  NF  propellants  may  pro¬ 
duce  a  100%  gaseous  exhaust.  In  the  consideration  of  the  N? 
propellants,  we  must  balance  the  potentially  poor  temperature 
sensitivity  against  the  highly  desirable  exhaust  which  is  free 
of  solid  particles. 

The  marked  difference  in  the  *-emperature  sensitivity 
between  the  NF  and  the  boron  compound  propellants  can  be  traced 
to  the  different  mechanisms  of  their  burning.  The  NF*s  have 
a  high  heat  release  since  they  burn  as  energetic  binders.  The 
binder  in  boron  compound  propellants  is  pyrolyzed  as  a  conven¬ 
tional  hydrocarbon  binder  and  therefore  reacts  in  the  gas 
phase.  It  is  expected  that  the  boron  compound  propellants 
have  surface  temperature  similar  to  the  hydrocarbon  binders 
and  that  the  NF  propellants  would  have  a  lower  burning  surface 
temperature  corresponding  to  that  of  dcuble-base  propellants. 

These  trendy  are  consistent  with  the  trends  predicted  by  the 
KTSS  model2  f  which  showed  that  increasing  the  heat  release 
or.  the  surface  <e,g.f  using  an  energetic  binder!  and  lowering 
the  surface  temperature  tend  to  increase  the  temperature 
sensitivity  of  burning  rate. 

*  Since  the  NF  propellants  have  been  found  to  be  particle  rree  in 
the  Hit  operating  pressure  range,  it  muy  be  necessary  to  include 
an  additive  such  as  carbon  black  as  a  suppressant  for  combus¬ 
tion  instability. 
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Another  well  characterized  family  of  propellants,  the 
composite  modified  double-base  propellants,  can  also  meet  the 
Hit  burning  rate  requirements  by  using  moderately  small  A?. 

The  conventional  ammonium  perchlorate  composite  propel¬ 
lants  can  meet  the  burning  rate  requirements  only  if  they  are 
catalyzed  with  metal  oxide  producing  burning  rate  catalyst. 
However,  their  high  specific  impulses  depend  on  high  percent¬ 
ages  of  ammonium  perchlorate  which  result  in  very  viscous 
propellants  that  are  difficult  to  cast.  There  is  an  interest¬ 
ing  trade-off  bee  *een  the  viscosity  of  high  AP  percentage 
propellants  and  the  easier  to  cast  propellant  grains  that  take 
advantage  of  the  higher  burning  rate  propellants.  When  iron 
compounds  are  added  as  burning  rate  catalyst,  metal  oxides 
will  form  in  the  exhaust.  For  example,  to  relate  the  percent¬ 
age  of  ferrocene  added  to  a  propellant  to  the  percentage  of 
metal  oxides  xn  the  exhaust,  an  approximate  calculation 
indicates  that  2%  ferrocene  will  produce  0.4%  Fe  O  in  the 
exhaust.  2  3 

Examining  the  various  propellant  systems  in  terms  of 
storability,  a  large  amount  of  data  exists  on  the  primary 
propellant  ingredients  such  as  ammonium  perchlorate  and 
binders.  However,  several  of  the  ingredients  such  as 
nitroglycerine  and  ferrocene  plasticizers  must  be  carefully 
examined  before  tnev  can  be  considered  for  either  vacuum 
storage  or  storage  in  fiber  glass  type  motor  cases. 


MOTOR  DESIGN  AND  TRADE-OFF  CONSIDERATIONS 


Baseline  Motor  Dtcsiqn 

The  design  of  a  baseline  no  tor  is  a  key  step  in  defining 
the  problems  thac  require  further  study.  The  propellant 
requirements  cannot  be  established  until  trade-offs  between 
items  such  as  burning  rate,  chamber  pressure  specif ic 

impulse,  and  case  diameter  are  evaluated  ana  understood*  The 
information  for  Table  3  was  extracted  from  Ref.14  and  from 
conversations  with  ABMDA  and  LTV  personnel.  Table  3  is  a  use¬ 
ful  starting  point  for  these  studies,  however,  the  relative 
importance  of  the  various  items  in  the  specification  must  be 
established  before  a  thorough  trade-off  analysis  can  be 
completed.  For  example,  it  may  b-.  possible  to  reduce  action 
time  by  reducing  the  mass  fraction.  Also,  the  mass  fraction 
cs.n  be  increased  by  using  a  noz2le  with  an  elxptical  exit  but 
only  at  an  added  fabrication  expense  to  the  end-item  assembly  - 

Figure  17  is  a  schematic  drawing  which  shows  the  case  and 
nozzle  configuration  which  was  considered.  The  equations  for 
the  dinvensions ,  weights,  and  volumes  of  the  case  and  nozzle 
components  have  been  programmed  for  computer  solutions.  The 
solutions  are  formulated  so  that  the  primary  constraints  are 
case  diameter,  port  to  throat  area,  and  total  impulse.  Light¬ 
weight  nozzles  are  considered.  For  example.  Figure  18  shows 
a  submerged  type  of  nozzle  that  has  a  throat  insert  and 
refractory  insulation  on  the  exit  cone.  No  attempt  is  made 
to  include  all  of  the  details  of  items  such  as  flanges  and 
welds . 

The  range  of  designs  which  can  be  considered  is  constrained 
by  three  special  features  of  the  Kit  impulsive  thruster:  the 
very  short  action  time,  the  high  length  to  diameter  r~tio,  and 
the  center  ’ented  chamber.  A  range  of  workable  motor  designs 
were  arrived  at  by  an  iterative  procedure.  We  will  discuss 
the  calculations  and  results  that  were  used  to  make  our 
decisions  by  making  comparisons  to  the  baseline  motor  described 
in  Table  4.  The  results  of  Fig. 19  relate  the  range  of  burning 
rates  required  to  meet  the  action  time  requirement  to  the 
available  web  fractions  shown  on  Fig. 20  The  summary  on  Fig. 20 
indicates  that  star  type  of  propellant  configuration  can  be 
used  with  a  case  diameter  of  0.30  inches  and  a  burning  rate 
of  10.0  in/sec.  {The  bases  for  selecting  a  case  diameter  of 
0.30  inches  will  be  discussed  when  Fig. 23  is  discussed.)  If  a 
higher  burning  rate  propellant  is  selected,  e.g..  13  in/sec. 
a  wago’.whael  type  of  propellant  configuration  will  be  required. 
As  indicated  on  Fig. 31  wago.uvheel  configurations  generally 
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TAB  LF  1 

EXAMPLE  SPECIFICATION  FOR  MICRO-ROCKET 

Units  are  inches,  seconds,  pounds  mass,  pounds  force,  and 
degrees  P . 


Total  impulse,  IT,  Ibf  sec 

Action  time,  fc 

Peak  thrust,  Ibf 

Thrust  misalignment,  3  o,  deg 

Outside  diameter  of  case,  U  ,  in 

co 

Minimum  design  factor,  u 

Maximum  temperature  on  outside 
of  case,  F 

Mass  fraction  of  complete  motor 

Overall  length,  m 

Ignition  delay  to  first  thrust ,  sec 

Variation  on  above  ignition 
delay,  3o,  5 

Variation  of  centroid  of 
pulse,  3c?,  sec 

Metal  oxides  and  other  solids  in 
exhaust,  % 

Operating  conditions 

Motor  temperature-  F 

Axial  acceleration  g's 

Spin  rate  about  central  axis  rps 


Baseline 

3.S 

0 . 005-0 . 008 
700.0 
0.25 

0.25-0.45 

1.3 

200.0 

2:0.49 

>6.0 

<0.0015 

<+20 

<0.00125 

<0.1 

vacuum 

50-70 

120 

25 
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Submot'gcid  Nos!*l«  Configuration  and  Dluvanwlom  Conaidorod  by 
Program  (Di monnionw  Aro  Not  .1  ?t  Proportion  for  Micromotor 
Configuration) . 
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_ _ ScnCTTQSS 

AMONG  ^G?=;LIAjfT  PAPAHTTERS 
{Baseline  Design  i?e.  1) 


units  are  inches,  seconds,  pounds  sass,  pounds  force,  pounds 
per  square  inch,  and  degree  F. 

Overall  Per forr-a nee  Paxasetcrs 

*  Total  impulse  during  t.  3.59 

Tise,  web  burning,  i  .  0.0069 

* ,  action,  t  - 

' ,  total  burning 


*  Pressure,  average  during,  t  10,000. 

,  saxissaa  instantas sous  ~I2,OOC. 

♦  Thrust,  average  during  t  ,  U>t  522. 

Overall  design  ,'acror,  y  1.4 

Total  irpulse/sotoi  «?:ght,  B2.4 

Fuel  weight/saofccr  weight,  0.345 


Propellant 

Type 

Xonsataliized,  A?  cos- 
posit  e  with  energetic 
binder. 

Specific  impulse. reference 

243.  1 

"  ",  delivered,  ISJ>D 

250.63 

Characteristic  velocity,  ft/sec 

5000. 

Suraing  rate  at  10,000  psi,in/sec 

10.0 

Burning  rate  exponent 

0.55 

temperature  sensitivity  of  burning  rate  0.0015 

at  constant  pressure,  -/? 


temperature  sensitivity  of  burning  rate  0.0034 

at  constant  A ^/A,. 


■* 

A  prescribed  value 
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TABLE 


•Ccr.t .  } 


Overall  Dissensions 
Case  length 

*  Case  outside  diameter 

Length  of  exit  cone  beyond  case 
Length  to  di ace ter  ratio 
Overall  Weights 
Propellant 
Case 

Case  liner  and  insulation 

Scrzle 

Igniter 

Motor 


5.798 
0.30 
0.  5S 
19.3:1 


0.015( 


0.020? 


<  0.001 


0.0076 


0.0436 


Propellant  Geometrical  Parameters 


Configuration 


star 


naming 


Outside  dianeter  of  propellant 
Characteristic  length,  initial,  L- 
Cross-sectional  loading  density 
Sliver  fraction 
Web  thickness 

Maxisus  Hach  nutber  in  port  after  norrle 
opens 

Liner  thickness 

Effective  port  area  to  threat  ares 
after  10%  of  veb  has  burned 

Initial  port  area  to  throat  area 


0.276 

3.072 

0.716 

0.05 

0.069 


<0.6 


~  0.002 


1.23 


l.GS 


—  *5 


?iSI-E  4  (Cent 


Chsaoet 

Type 

Material 

,  strength 
SesBiasl  thscfaass 
Insniatic-ri 


level ,  psia 


Hollcv  tiibisc  vith  b< 

635  cap 

Haragin:  stasis  Grace  2: 
240, GOD 
0.010 

Case  wall  arc  end  cans 


Case  wall  thickness  {0.010  is 
sdniszst  valce  frees  sats- 
csssidsratios) 


0.010 


Settle 


Type 

rsv'ir'*  exit  wSdr 

refractory  coating 

Sx^rsics  Ratio 

3.74 

Divergence  angle 

15.0 

Sxit  Dissaeter 

0.60 

Throat  diaseter 

0.203 

Shell  .Material 

Kara ginc  steel.  Grade  2>0 

'  strength 

level 

4*rV  t  vVt- 

Insert  material 

Refractory  coating 

Shell  wall  thickness 

0.010 

m  *  * 

at  exit 

0.005 

6 
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8  10 

BURNING  RATE,  r,  i^sec 

feet  of  burning  rate  and  chamber  pressure 
.  burning  surface  area  and'  web  fraction  . 

'  baseline  motor. 


DENDRITE 

(FORKED 

WAGONWHEEL) 


WAGONrfxiEEL 


STAR 


SLOTTED 

TUBE 


CONFIGURATION  WAGONWHEEL  FORKED  STAR,  STAR, 

_ _ _ WAGONWHEEL  IDEAL  THICK 

WEa/RADIorB  0-12-0.17  0.20-0.25  0.3-0~4  5-0.8  lo.S-0.' 


GROSS  VOLUME  LOADI1 

SLIVER  (IDEAL) 

CONFIGURATION 

EFFICIENCY^ 

MAXIMUM  L/D  FOR 
NEUTRALITY* 

CASE  EXPOSURE 


75%  85%  95% 


93%  95%  98% 


UNLIMITED  5: 1-8: 1  4:1-5:1  L 1-4:1  b:l-3; 


NONE  TO  WEI  NONE  TO  WEI  NONE  TO  NONE  TO  PRIOR 

WEB  WEB  TO  WEE 


COMMENTS  PERTAINING  LOW  VOLUME  GOOD 
70  HIT  LOADING 


GOOD  LOW  LOW 
L/D  l/d 


*L/D  IS  BASED  ON  U./2  OF  MOTOR!  LENGTH 


Fig.  20  Relative  web  *pe ctrurc  and  attributes  of  neutral  hurHr, 
grain  design*  (after  BiUheimer)  *al  bu‘ni^ 


CONFIGURATION  EFFICIENCY 


52 


Fig.  21  Configuration  efficienty  as  a  function  of  volumetric 
loading  factor  for  grain  designs  of  various  typical 
web  thicknesses  (after  Billheiraer) 


-  D  2  ~ 


have  lower  volumetric  loading  densities  and  configuration 
efficiencies  than  starpoint  configurations.35  Thus  using 
the  lower  burning  rate  of  10.0  in/sec  does  not  penalize  the 
system  and  gives  us  more  latitude  in  selecting  the  propellant 
formulation.  Another  consideration  for  a  neutral  burning 
propellant  configuration  is  the  burning  surface  perimeter. 

Figure  22  shows  the  range  of  the  ratio  of  burning  perimeter 
to  the  circumference  of  the  grain.  The  baseline  design  is 
in  the  proper  range  to  achieve  a  neutral  burning  propellant 
configuration.  However,  as  the  buring  rate  increases,  it 
becomes  progressively  more  difficult  to  achieve  a  neutral 
burning  propellant  configuration.  (Neutral  burning  means 
that  the  burning  surface  area  stays  reasonably  constant  as 
the  propellant  burns.) 

One  of  the  constraints  on  the  system  is  that  the  motor 
must  exceed  a  minimum  length.  Figure  23  shows  how  a  range 
of  lengths  affects  volumetric  loading  density.  If  the  case 
diameter  is  too  small,  e.g.,  0.25  in.,  the  loading  density 
is  very  small  since  the  cross-sectional  area  of  the  pert  must 
be  larger  than  the  nozzle  throat  area  (e.g.,  A  ort/At>l .0 ) . 

The  rapidly  increasing  Mach  numbers  of  the  internal  flows  as 
A  t/At  approaches  unity  produce  serious  performance  losses 
such  as  pressure  drops  in  the  direction  of  flow,  nonuniform 
burnout  of  the  propellant  web,  and  thrust  misalignment.  The 
larger  diameter  motors  are  more  inefficient  since  volumetric 
loading  densities  are  too  low. 

Even  if  a  relatively  small  initial  port  area  to  throat 
area  ratio  is  used,  the  high  degree  of  surface  blowing  pro¬ 
duced  by  the  very  high  burning  rate  propellants  eliminate 
serious  erosive  burning  effects.  This  observation  is  supported 
by  the  recent  calculations  of ^illoughby ‘  and  was  reviewed  in 
terms  of  Saderholm's  criterion 

Even  though  the  action  time  is  short,  the  very  high  opera¬ 
ting  pressures  being  considered  complicate  the  nozzle  design. 
Several  series  of  heat  transfer  calculations  were  conducted 
to  evaluate  the  nozzle  insulation  requirements.  Two  design 
situations  were  considered: 

1)  an  uninsulated  steel  throat  section  and  exit  cone 

2)  a  0.010  in  of  A^Oq  type  of  refractory  coating  on 
a  steel  throat  section. 

The  convective  heat  transfer  coefficients  were  estimated 
using  the  equations  of  Ref.  36,  and  the  baseline  propellant 
gas  properties  Figure  24  shows  that  situation  1  is  marginal 
for  chamber  pressures  10,000  psi  and  unsatisfactory  for  15,000 
psi.  Situation  2  appears  to  be  generally  satisfactory  since 
the  surface  temperature  is  well  below  the  melting  point  of 
the  refractory  coating  and  since  there  is  only  a  slight  increase 
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pig.  23 


ca“®  aiamete*  and  volumetric 

0.5) 

£* 


400  800  1200  1600 


M  TEMPERATURE  OF 
SEVERE  EROSION 


24  Surface  temperature  of  nozzle  throat: 
’uninsulated  steel  and  a  refractory 
coating  on  a  steel  shell. 
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in  temperature  at  the  refractory  coatir.g/steel  interface. 
Calculations  were  also  carried  out  for  the  exit  cone  (Fig.  25). 

The  ratio  .... /W’  is  a  good  performance  indicator 
for  the  motors  w4'are  studying.  The  ratio  W  .  /Wwo 
(used  in  Ref.  34)  is  not  as  revealing  since  ?£°§oes^not  fully 
account  for  the  importance  of  specific  impulse.  Figure  26 
shows  that  the  baseline  case  diameter  does  not  produce  the 
most  efficient  motor.  However,  the  L/D  requirement  limits  us 
to  a  diameter  of  approximately  0.30  inches.  As  shown  in  Fig.  26, 
the  nozzle  weight  is  a  large  fraction  of  the  overall  motor 
weight.  Because  the  chamber  pressure  is  high  and  the  ambient 
conditions  are  vacuum,  large  increases  (e.g.,  30  sec.)  in  Isc<j 
can  be  achieved  by  increasing  the  exit  diameter.  To  limit 
the  vibrations  in  the  system,  it  is  necessary  to  use  a  rather  ^ 
stiff  nozzle  ring  and  to  firmly  secure  the  exit  cones  in  place3  . 
For  these  reasons  the  exit  cone  shell  thickness  (~R4  or.  Fig.  18) 
was  maintained  at  a  minimum  value  of  0.005  inches.  ‘  The  net 
result  is  that  any  performance  increase  obtained  by  a  larger 
exit  cone  is  nullified  by  the  increased  weight  of  the  exit 
cone.  This  trend  is  shown  in  Fig.  26  when  the  exit  diameter 
is  decreased  from  1.5  to  1.0  inches.  Figure  27  shows  that  a 
good  value  of  DE  is  approximately  0.6.  Another  important 
result  shown  on  Fig.  27  is  the  effect  of  increasing  chamber 
pressure.  The  increased  ISDCj  and  loading  density  correspond¬ 
ing  to  the  higher  chamber  pressure  are  nullified  by  the 
increase  in  case  weight. 

The  interpretation  of  results  of  Figs,  19  through  27 
is  a  first  step  in  defining  the  operating  range  and  propel¬ 
lants.  However,  interpeting  the  results  has  to  be  evaluated 
in  terms  of  items  such  as:  the  basis  for  selecting  a  particular 
motor  length,  the  relative  importance  of  It/wmo  terms  of 
material  costs  and  ease  of  handling,  the  penalty  for  increas¬ 
ing  action  time,  and  conditions  under  which  increasing  action 
time  while  increasing  I«i>/WWq  is  a  net  gain. 

Calculated  Pressure  responses 

The  analysis  section  emphasized  the  necessity  of  accurately 
treating  the  transient  burnina  rate  and  the  interactions  between 
combustion  and  gas  dynamics.  The  results *in  Figs.  28  and  29 
are  typical  of  the  results  which  are  obtained  from  the  tran- 
ient  gas  dynamics  and  burning  rate  model.  The  uniform  p  vs  t 
during  the  first  two  milliseconds  (shown  on  Fig.  29)  was 
obtained  after  conducting  several  computer  experiments  where 
igniter  mass  flux,  igniter  heat  flux,  nozzle  closure  blow-out 
pressure,  and  propellants  were  varied.  Figure  28  shows  some 
of  the  details  considered  in  selecting  propellants,  operating 
ranges,  and  igniters.  Note  that  at  approximately  0.5  milli¬ 
seconds,  the  propellant  becomes  fully  ignited  and  any  addi¬ 
tional  contribution  from  the  igniter  is  not  required.  One 
item  that  is  necessary  to  achieve  precision  of  the  thrust 
versus  time  trace  is  to  find  the  domain  where  the  propellant 

*  These  results  are  for  the  design  described  in  Table  5  which 
will  be  discussed  in  the  section  entitled  "Performance  Analyzed 
in  Terms  of  F  vs  t . " 
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Fig. 25  Surface  temperature  of  nozzle  exit} 
uninsulated  steel  and 0.002  inches  of 
a  refractory  coating  on  a  steel  sheel 


UPPER  MAP  Dg=  1,0 
LOWER  IftP  *  1.5 


-*  BASELINE  MOTOR  EXCEPT  FOR  D£ 

Pig, 26  Parametric  map  of  effect  of  loading 
density,  case  diameter,  and  nozzle 
exit  diameter  on  motor  specific  impulse 


-*  BASZLIKE  MOfTOH. 


MJCRD-**DT03  ?£RTORKASC£  PARAKZzSRS 
{Basel!.’;?  Design  So.  2} 

Snits  are  inches ,  seconds ,  poonds  sass,  pounds  force,  pounds 
per  square  inch,  and  decree  ?. 


Overall  Performance  Par  acre -^ers 

Directed  total  isspa is  e/t  o  ta  1  impulse  ,  I  /i„„ 
?ise,  web  burning, 

*  ,  action,  t 

*  ,  total  burning  and  tail— off  to 
0.-5  05  i 

'CS 

*  Pressure,  jyssinal  during,  t 

,  aauci sssst  instantaneous 

*  Thrust ,  average  during  t  ,  ibf 

3 

Overall  design  factor. 

Total  inpalse/notcr  weight, 

Props  1  ia nt  weight /no  tor  weight,  Vilj 


3.32/4.0! 

C.005S 

0.0075 

0.014 
15,000 
<  10,000. 
525  . 

1.3 

32.3 
U.40 


Praoeilant 

Type  (So.  1  in  Table  3}  Sense  tall!  zed,  A?  ccssxos 

f  Specific  inpalse.  reference  delivered  228-7 

*  ",  delivered  1  . 

Spd 

Characteristic  velocity,  ft/sec 
Density, 

Burning  rate  at  10.OD0  psi.  in/sec 
3cm  me  rate  ea 


ou.u 


f.5»C 


•Peeper  a  cure  sensitivity  at  corn  mg 
at  constant  pressure,  r'^- 

Tesperature  sensitivity  of  burning 
at  constant  A. /A_ 


rate 


A  prescribed  value 

Per iim.cce  not  degraded  for  losses  resulting  r'r.o  center 
vented  nozzle 
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TABLE  5  (Cont.) 


Overall  Dimensions 

Case  length 

7.91 

*  Case  outside  diameter 

0.31 

Length  of  exit  cone  beyond  case 

0.53 

Length  to  diameter  ratio 

25.5:1 

Overall  Weights 

Propellant,  w 

prop 

Case 

0.0160 

0.0217 

Case  liner  and  insulation 

None 

Nozzle 

0.0038 

Igniter 

- 

MOt°r'  WMO 

Propellant  Geometrical  Parameters 

0.0415 

Configuration 

Internal  burning 
double  web 

Outside  diameter  of  propellant 

0.300 

Characteristic  length,  initial,  L* 

11.22 

Cross-sectional  loading  density 

0.53 

Sliver  fraction 

0.05 

Web  thickness 

0.073 

Maximum  Mach  number  in  port  after  nozzle 
opens 

0.21 

Effective  port  area  to  throat  area  after 
10%  of  web  has  burned 

1.5 

Initial  port  area  to  throat  area 

1.4 

'ffflRJT  IP* 
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TABLE  5  (  cont.) 


Chamber 

Type 

Material 

'*  ,  strength  level,,  psia 
Nominal  thickness 
insulation 

Case  wall  thickness  (0.005  is 
minimum  value  from  manu¬ 
facturing  consideration) 


Hollow  cubing  with 
bonded  end  cap 

Maraging  steel 

300,000 

0.010 

Case  wall  not  insulated; 
end  caps  lined 

0.030 


Nozzle 

Type 

Expansion  Ratio 
Divergence  angle 

°*  eliptical 

Ihroat  diameter 
Shell  Material 

"  ,  strength  level,  psi 
Insert  material 

Shell  wall  thickness  at  insert 


Type 


Propellant  weight 
Duration 

Heat  flux  to  surface 
Propellant  type 


Submerged;  conical  exit 
cone;  refractory  coating 

12.5 

15.0 

0.60 

0.168 

Maraging  steel.  Grade  250 
240,000 

Refractory  coating 

0.014 

0.005 

Squib  at  one  end 
0. 10  g 
0 . 001  sec 
200  cai/cn^sec 


No,  1  in  Table  2 


-  OH 


RESULTS  ARE  FOR  MOTOR  SIMILAR  TO  MOTOR 
DESCRIBED  IN  TABLE  5  EXCEPT  NOMINAL 
OPERATING  PRESSURE  IS  10,000  PSI  RATHER 


0.0  0.002  0.004  0.006  0.008 


TIME  FROM  ONSET  OF  IGNITION  STIMULUS,  sec 

Fig.  28  Chamber  temperature  and  burning  rate 
dynamics  during  motor  operation 


STATIC  CHAMWER  PRESSURE 


Fig'  29  Head-end  and  nozzle-end  pressures 


STATIC  PRESSURE/FINAL  EQUILIBRIUM  STATIC  PRESSURE 


RESULTS  ARE  FOR  MOTOR  SIMILAR  TO  MOTOR  DESCRIBED  IN 
TABLE  5  EXCEPT  NOMINAL  OPERATING  PRESSURE  IS  10 , 000 
PSI  RATHER  THAN  15,000. 


TIME,  t,  sec 


Fig.  30  Effect  of  decreasing  throat  diameter  on  pressure 
versus  time  history 
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becomes  fully  ignited  in  a  uniform  manner.  For  example,  if 
the  igniter  heat  flux  is  too  great  the  mass  contribution 
from  the  propellant  will  depend  on  the  rate  at  which  the 
igniter  volatilizes  the  propellant  surface  rather  than  on 
the  inherently  more  stable  process  of  self  sustaining  pro¬ 
pellant  combustion.  Relatively  small  igniters  produce 
extremely  rapid  flame  spreading  rates  along  the  propellant 
surface  since  the  conditions  of  high  burning  rate  propellants, 
high  pressure,  and  low  A _/At  combine  to  increase  the  heat 
transfer  to  the  propellant 


-3  T-l  9 


The  decrease  between  the  head-end  and  the  nozzle-end 
pressure  is  a  result  of  the  acceleration  of  the  gases  along 
the  port.  It  may  be  desirable  to  reduce  the  pressure  drop 
in  order  to  achieve  a  more  efficient  motor  case  design. 

Figure  30  shows  that  decreasing  the  nozzle  throat  diameter 
so  that  the  motor  operates  at  15,000  psi  rather  than  10,000 
psi  increases  A^rt/At  and  thus  decreases  the  pressure  drop 
along  the  port.*in  the  15,000  psi  case,  the  nozzle  closure  was 
set  to  open  at  95%  of  the  operating  pressure  and  in  the  10,000 
psi  case  the  nozzle  closure  was  set  for  75%  of  the  operating 
pressure. 

While  it  was  not  possible  in  this  study  to  evaluate  all 
of  the  trade-offs  in  the  motor  selection  process,  several  of 
the  design  parameters  used  to  achieve  the  wDrQD/WMQ  *-n  the 
specification  (Table  3)  indicate  clearly  thSt°TabXe  3  is 
presently  unrealistic.  Table  5  describes  a  motor  that  satis¬ 
fies  many  of  the  requirements  in  the  specification.  However, 
the  basis  for  several  of  the  parameters  is  questionable.  To 
achieve  the  lew  inert  hardware  weight  miniature  seals  and 
flanges  were  assumed  and  very  high  strength  steel  was  used. 

A  higher  propellant  loading  density  could  be  obtained  by  using 
very  low  port  area  to  threat  area  ratio  which  results  in 
high  Mach  numbers  (0.3)  in  the  port.  This  is  a  particularly 
serious  situation  since  the  two  opposing  high  Mach  number 
streams  must  be  turned  90  degrees.  The  decrease  in  specific 
impulse  resulting  from  turning  the  two  opposing  flows  and 
the  short  flow  path  after  the  turn  has  net  been  established - 
To  proceed  into  the  next  phase  of  a  motor  development  program 
without  more  information  on  the  effects  of  opposing  flows 
could  result  in  the  selection  of  motors  that  are  light  weight 
but  with  inherently  unsatisfactory  thrust  alignment  and  non- 
repeatabxe  thrust  versus  time  traces.  In  other  words, 
adequate  information  must  be  obtained  for  establishing  at 
what  point  decreases  in  motor  weight  actually  degrade  the 
overall  s/stem  performance. 


EXHAUST  PLUME  CONSIDERATIONS 


The  exhaust  plume  considerations  are  complex  and  involve 
several  overlapping  objectives  and  requirements  which  will 
occur  during  the  development,  demonstration,  and  operational 
phases.  The  literature  on  contamination  from  exhaust  impinge- 
ment,,“<*9  and  flow  fields  is  only  partially  applicable. 

Our  interpretation  of  a  reasonable  course  of  action  is  to  use 
as  a  goal  the  requirements  as  set  down  in  the  specification 
(Table  3)  but  not  to  penalize  the  performance  of  prototype 
impulsive  thrusters  by  items  that  are  not  important  to  the 
prototype  tests.  In  this  manner,  the  system  can  be  developed 
and  tested  while  allowing  for  improvements  and  changes  that 
will  ^ccur  during  the  development  and  demonstration  phases. 

For  example,  during  the  development  and  demonstration  phases 
long  term  storage  is  not  required;  however,  reproducibility 
of  the  centroid  of  thrust  versus  time  traces  is  very  important. 
Similarly,  many  of  the  tests  can  be  conducted  using  propellants 
that  produce  submicron  size,  small  amounts  of  metal  oxide 
particles  (<  2%)  in  the  exhaust.  Since  the  detrimental  effects 
of  metal  oxides  have  not  been  established,  studies  to  evaluate 
metal  oxide  effects  can  be  conducted  concurrently  with  test 
programs  that  use  prototype  versions  of  the  impulse  thruster. 

For  example,  one  of  the  items  to  be  evaluated  is  the  size  of 
the  metal  oxide  particles  (in  the  exhaust  plume)  produced  by 
iron  compound ,  burning  rate  catalysts.  If  the  iron  oxide  does 
not  agglomerate  into  large  particles,  then  the  thermal  response 
of  the  iron  oxide  may  closely  follow  the  thermal  response  of 
the  expanding  exhaust  gases.  It  appears  likely  that  the  present 
Hit  impulsive  thruster  requirements  can  be  met  with  a  propel¬ 
lant  that  does  not  produce  metal  oxides.  However,  improved 
motor  performance  can  be  achieved  if  the  present  limitation 
on  iron  compound /burning  rate  catalysts  can  be  relaxed. 

Table  6  lists  some  of  the  developmental  test  conditions 
that  apply  to  metal  oxides  in  the  exhaust.  It  is  important 
to  note  that  the  three  test  conditions  are  very  different  from 
each  other.  (We  realize  that  Table  5  deals  with  only  one 
aspect  of  the  exhaust  plume  considerations.)  In  all  situa¬ 
tions,  it  is  unlikely  that  appreciable  quantities  of  the 
relatively  heavy  metal  oxide  particles  will  come  into  contact 
with  the  system.  A  potential  source  for  concern  is  the  accu¬ 
mulation  of  gaseous  exhaust  products  during  the  tests  under 
vacuum  conditions.  In  that  case,  the  lighter  molecules  may 
be  the  most  likely  to  come  into  contact  with  system. 

Exhaust  Emissions 


The  reasons  for  eliminating  metal  oxides  from  the  exhaust 
can  stem  from  several  requirements.  One  requirement  is  that 
as  the  hot  metal  oxides  are  exhausted,  they  are  slow  to  cool 
dewn  under  vacuum  conditions  and  therefore  can  be  detected  by 
various  heat  sensors.  Internal  ballistic  studies  indicate  that 
the  exhaust  pressure  will  be  greater  than  50  psi  and  thus  the 
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TABLE  6 


COMMENTS  OK  DEVELOPMENTAL  TEST  CONDITIONS  THAT  RELATE  TO 
METAL  OX IDES  IN  THE  EXHAUST 


Tests  at  normal  Tests  under  j  Flight 

.  ambient  conditions  !  vacuum  conditions  •!  tests 

- p - 1 - S - 

Environment  *  Ig  ]  la  I 

,  1  atmosphere  \  -‘0.05  atmosphere  i  ^0.05  atmo- 

60  to  80  F  :  !  sphere 


Gaseous  exhaust 
considerations 


Smoke  aay  be  removed  Exhaust  gases  can 
by  low  velocity  air  accumulate  in 

flowing  past  motor.  vicinity  of 

system. 


System  is 
accelerated 
away  from 
exhaust 
gases . 


Metal  oxide 
considerations 


For  a  motor  firec 
downward  gravity  and 
momentum  of  exhaust 
products  act  so  as 
to  prevent  metal 
oxides  from  accumu¬ 
lating  in  the  vicin¬ 
ity  of  the  system. 


For  a  motor  fired 
downward  gravity 
and  momentum  of 
exhaust  products 
act  so  as  to  pre¬ 
vent  metal  oxides 
from  accumulating 
in  the  vicinity 
of  the  system. 
However,  some  of 
the  lighter  gases 
may  diffuse  to¬ 
ward  the  system. 


I 


I 

f 


1 

I 

j 


The  individ¬ 
ual  metal 
oxide  parti¬ 
cles  in  the 
exhaust  have 
much  greater 
momentum 
than  the 
molecules  of 
the  gaseous 
products  and 
thus  are  no*- 
likely  to  be 
entrained  by 
the  system. 


Storage  period 


1  yr. 


Operating  50  to  8C  F 

temperature 


-  1  yr. 


•1  yr. 


Not  established  Not 

I  established 


plume  will  be  relatively  dense.  An  all  gaseous  exr.aust  will 
rapidly  dissipate  and  leave  a  very  low  thermal  signature. 

Radar  return  can  i  •»  greatly  influenced  by  the  ffsa 
electrons  in  the  exhaust.  NF  propellants  produce  C},  to 
attract  the  free  electrons.  The  CIO  in  the  exhaust  is  not 
disturbed  by  the  radar  because  of  its  relatively  high  mass. 
However,  the  low  mass  free  electrons  in  the  exhaust  are 
excited  by  the  radar  and  will  reflect  radar  waves  and  diminish 
the  transmission  through  the  exhaust  products.  The  serious¬ 
ness  of  this  depends  on  whether  the  radar  is  communicating 
with  the  vehicle  or  whether  the  radar  is  trying  to  locate  the 
vehicle.  Propellants  containin-  aluminum  will  have  a  higher 
flame  temperature  which  vill  pi^^uce  more  ionization,  and 
thus  more  free  electrons  and  higher  thermal  emissions. 

The  wave  lengths  of  radar  transmissions  are  large 
compared  to  the  micron  size  of  the  exhaust  products.  Accord¬ 
ingly,  the  scattering  cross  sections  of  the  particles  are  much 
smaller  than  the  actual  cross  sections.  Thus,  the  Hit  system 
operates  in  a  good  regime  as  far  as  radar  detection  is  concerned. 
However,  it  is  not  in  a  good  regime  as  ;ar  as  near  infra-red  is 
concerned  because  the  particle  diameter  is  an  appreciable 
fraction  of  (or  larger  than)  the  wave  length  of  the  radiation. 

The  benefit  of  replacing  A?  with  KMX  in  smokeless  pro¬ 
pellants  at  sea  level  conditions  is  that  the  concentration  of 
HC1  is  reduced.  KC1  promotes  the  condensation  of  water  by 
reducing  the  vapor  pressure  of  small  droplets.  Again,  because 
of  vacuum  conditions,  the  likelihood  of  condensation  of  water 
vapor  is  remote  ever,  with  the  oresence  of  a  high  concentration 
of  HCI . 

Back  Flow  of  Exhaust  Gases 


As  the  under-expanded  gases  flow  past  the  exit  plane  of 
the  norzle ,  they  undergo  further  exoansion  and  turn  away  from 
the  centerline  of  the  nozzle.  This  will  be  referred  to  as  thr 
plume  interference  problem.  Under  the  conditions  of  sufficiently 
low  back  pressure  and  high  exit  pressure,  it  is  possible  (in 
principle)  for  part  of  the  exhaust  gases  to  turn  100  degrees 
and  impinge  on  the  motor  case.  This  impingement  (also  referred 
to  as  spillage)  can  have  complicating  effects  and  should  be 
avoided . 


The  turning  angle  of  the  exhaust  gases  flowing  past,  the 
nozzle  exit  plane  {i.e.,  shown  in  Fig.  31)  can  be  approxi¬ 
mated  in  terms  of  the  Prandfcl-Meyer  angle.  This  turning  angle 


is  a 


irst  approximation  tc  the  limiting  stream  line  of  the 
exhaust  plume.  As  a  general  requirement,  should  be  less 
than  90s.  Figure  32  illustrates  those  items  which  tend  to 
reduce  &*•  Increasing  the  diameter  of  the  exit  plane  (i.e., 
the  expansion  ratio)  and  the  ratio  of  specific  heats  of  the 


exhaust  gases  decrease  the  turning  angle.  The  effective  cross- 
sectional  area  of  the  exit  plane  can  Be  increased  (within  the 

present  envelope  constraints)  by  using  an  exit  cone  with  an 
eliptical  cross-section.  Indeed,  losses  resulting  from  the 
complex  flow  patterns  in  the  nozzle  resulting  from  the  center- 
vented  configuration,  may  be  reduced  by  using  a  throat  section 
with  an  eliptical  cross-section.  Such  a  nozzle  should  have  a 
more  uniform  pressure  distribution  at  the  exit  plane.  Also,-  as 
indicated  on  Fig.  31,  reducing  the  nozzle  exit  divergence  angle 
has  the  effect  of  straightening  the  flow.  It  is  important  to 
note  that  under  hard  vacuum  conditions  vehicle  movement  during 
motor  operation  has  no  appreciable  effect  on  the  shape  of  the 
exhaust  plume. 

A  series  of  calculations  was  performed  to  investigate 
those  items  which  tend  to  zeduca  the  likelihood  of  8^  being 
greater  than  90°.  As  shown  in  Fig.  33a  increasing  the  nozzle 
exit  diameter  decreases  IT  V;efc/WM''  ^i^ce  the  increase  in 
delivered  specific  impulse' iFig.  *^3b)  does  not  off-set  the 
increase  in  nozzle  weight  Fig.  33c.  As  indicated  in  Fig.  32, 
the  exhaust  gas  composition  has  a  large  influence  on  8^.  For 
very  short  nozzles,  the  question  naturally  arises  as  to 
whether  there  is  sufficient  time  for  shifting  composition  to 
be  valid.  Figure  33d  indicates  that  the  extent  to  which  the 
composition  does  not  shift  as  the  gas  flows  through  the  nozzle 
has  an  appreciable  effect  on  the  conditions  at  the  exit  plane. 

Figure  34  shows  the  combined  effects  of  tailoring  the 
propellant  so  that  it  has  a  higher  ratio  of  specific  heats 
and  increasing  the  expansion  ratio,  e.  Any  change  that  shifts 
the  operating  point  toward  the  upper  right  hand  corner  of  the 
figure  decreases  the  likelihood  of  spillage.  (The  theoretical 
boundary  of  no  spillage  is  shown  by  the  row  of  arrows) .  Since 
there  are  envelope  constraints  on  DE,e  can  be  increased  by 
decreasing  Afc  and  thereby  increasing  the  operating  pressure. 
Previous  discussions  pointed  out  that  if  an  AP  composite  pro¬ 
pellant  is  used,  the  tradeoffs  between  burning  rate, 

processability,  and  metal  oxides  in  the  exhaust  will  influence 
the  AP  percentage.  These  tradeoffs  are  further  complicated 
since  lowering  AP  percentage  significantly  reduces  the  like¬ 
lihood  of  spillage. 

The  foregoing  discussion  was  included  to  illustrate  the 
spillage  consideration  in  the  overall  motor  design  process. 
There  appears  to  be  reasonable  latitude  available  to  design 
around  the  occurrence  of  spillage. 
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Fig.  33  Motor  parameters  as  a  function  of  nozzle  exit  diameter 
(For  motor  described  by  Table  5) 
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PgRPORMANCE  ANALYZED  IN  TERMS  OF  F  vs  t 


The  transient  ballistic  performance  model  was  used  to 
study  how  various  motor  and  propellant  variables  influenced 
the  performance  of  the  motor  described  in  Table  5,  referred 
to  in  this  section  as  the  Baseline  Design.  The  results  are 
presented  by  comparing  thrust  versus  time  programs  (F  vs  t) . 
Since  the  Hit  impulsive  thrust  rotates  about  an  axis  parallel 
to  its  axial  center  line,  it  is  important  to  understand  how 
the  shape  of  F  vs  t  affects  the  directed  thrust  vector  (Tx) 
and  smear  factor  (It/It) •  These  terms  are  defined  and  dis¬ 
cussed  on  Page  34.  It  should  be  noted  that  the  calculated 
F  vs  t  results  are  idealized  compared  to  measured  F  vs  t  for 
systems  as  complex  as  the  Hit  impulsive  thruster.  The  re¬ 
sults  are  idealized  primarily  because  there  is  no  practical 
way  of  fully  accounting  for  the  nonuniform  burning  of  the 
propellant.  We  expect  that  the  sharply  defined  on-set  of 
tail-off  in  our  calculations  will  not  be  observed  even  in 
the  most  carefully  manufactured  Hit  impulsive  thrusters. 
However,  we  expect  that  the  extent  to  which  changes  in  motor 
and  propellant  variables  influence  both  measured  and  calcu¬ 
lated  F  vs  t  will  be  comparable. 

An  important  manufacturing  consideration  is  the  dimen¬ 
sional  tolerances  to  be  held  on  the  propellant  grain.  These 
tolerances  will  strongly  influence  the  final  selection  of 
the  type  of  grain  design  and  manufacturing  technique  to  be 
applied  (e.g.,  extrusion  of  an  internal  star  or  casting  of 
a  longitudinal  slot) .  As  an  example,  calculations  were  made 
to  simulate  the  effect  of  off-setting  a  star-point  mandrel 
by  O.OG35  in.  (about  5%  of  the  web  thickness) .  The  F  vs  t 
results  of  Fig.  35  show  that  the  effect  is  slight.  Even 
though  the  dimensions  cf  the  grains  are  small,  it  is  reason¬ 
able  to  expect  dimensional  consistency  to  within  ±  0.003  in. 

All  of  the  oerformance  nredictions  discussed  thus 
far  were  for  motors  with  nozzle  closures  designed  to  open 
very  near  to  the  motor  operating  pressure.  The  success  of 
nozzle  closures  (that  open  after  the  motor  is  near  its  oper¬ 
ating  pressure)  in  obtaining  reproducible  p  vs  t  histories 
has  been  conclusively  demonstrated  (e.g.,  the  micro-motor 
ballistic  test  device  developed  by  Rohm  and  Haas).  Since 
the  Hit  impulsive  thruster  is  to  be  stored  and  ignited  under 
vacuum  conditions,  some  type  of  nozzle  closure  is  essential. 
The  primary  questions  that  remain  are:  what  is  the  most 
desirable  opening  pressure  and  what  mechanical  design  is 
most  suitable  from  overall  system  considerations.  Figure 
36  shows  calculated  F  vs  t  for  nozzle  closures  with  removal 


NOZZLE  CLOSURE  OPENS 
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-  BASELINE  DESIGN,  TABLE  5 

_  SAME  AS  BASELINE  EXCEPT  SLIVER  IS  INCREASED  5% 

BY  OFF  SETTING  STAR -POINT  MANDREL  0.0035  INCHES 


Fig. 35 


Calculated  thrust  versus  time  programs  showing  how 
a  baseline  configuration  is  altered  by  core  misalignment 
which  increases  sliver  fraction  5 %. 
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-  BASELINE  DESIGN,:  TABLE  5 

-  SAME  AS  BASELINE  EXCEPT  PRESSURE  OF  NOZZLE  CLOSURE 

REMOVAL  IS  DECREASED  FROM  14.000  TO  5.000  PSI 

.  SAME  AS  BASELINE  EXCEPT  PRESSURE  OF  NOZZLE  CLOSURE 

r  REMOVAL  IS  DECREASED  FROM  14,000  PSI  TO  1.000  PSI 


0  0.002  0.004  0.006  0.008  0.010  0.012 

TIME  AFTER  ON  SET  OF  IGNITION  STIMULUS,  t,  SEC 

Fig.  36  Calculated  thrust  versus  ties  programs  showing  hew 
a  baseline  configuration  is  altered  by  decreasing 
pressure  of  noczle  closure  reaoval* 
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pressure  of  14,000  5,000,  and  1,000  psi.  As  the  nozzle 
closure  pressure  decreases, the  action  Lime  of  the  motor 
increases  and  the  F  vs  t  shape  becomes  progressively  more 
rounded.  On  first  appearances  one  might  expect  that  the 
more  rounded  F  vs  t  would  have  a  much  smaller  smear  factor 
(Tr/I?)  and  smaller  It.  However  since  the  motor  is  oper¬ 
ating  under  vacuum  conditions  with  a  prescribed  nozzle, 

I^pd  to  a  good  approximation  is  independent  of  pressure. 

Figure  3^  shows  that  approximately  2%  additional  losses  in 
Tr  coupled  with  thrust  misalignment  of  12  degrees  occurs  at 
the  specified  3c  limit  on  tcen.  Rounding  of  F  vs  t  as  in 
Fig.  36  results  m  only  slight  additional  losses:  an  addi¬ 
tional  0.5%  loss  in  directed  total  impulse  (Tr)  and  no  apprec¬ 
iable  change  in  the  degree  of  thrust  misalignment. 

Throughout  this  series  of  calculat ions,  it  was  noted 
consistently  that  small  departures  in  F  vs  t  had  insignifi¬ 
cant  effects  on  "tr.  This  observation  is  strong  evidence 
that  the  operating  point  of  the  Baseline  Design  has  been 
properly  selected.  As  an  example  of  the  effects  produced 
by  large  departures  ir.  the  Baseline  Deisng  F  vs  t,  F  vs  t 
was  calculated  for  a  nominal  pressure  of  10,000  psi  rather 
than  15,000  (see  Fig.  38).  This  change  increases  the  action 
time  by  0.0013  seconds  and  decreases  the  directed  total 
impulse  "fr  by  about  2%.* 

The  effect  of  decreasing  Ap0rt/At  has  been  discussed 
previously  in  terms  of  the  adverse  effects  that  the  higher 
internal  Mach  number  would  have  on  thrust  misalignment  and 
losses  produced  by  turning  the  port  flow  90  degrees.  Fig¬ 
ure  39  shows  a  F  vs  t  prediction  for  a  50%  reduction  in  the 
port  area.  Figure  40  shows  the  combined  effect  ot  decreasing 
Ap0rt  and  the  nozzle  closure  removc.1  pressure.  The  high  pres¬ 
sure  drop  along  the  grain  has  the  effect  of  rounding  off  F  vs  t 
and  the  increased  mass  addition  losses  decreased  the  total  im¬ 
pulse  by  3%.  The  turning  losses  were  not  accounted  for.  How¬ 
ever,  one  should  anticipate  that  turning  two  opposing  Mach 
0.46  streams  will  produce  serious  mixing  losses  at  the  con¬ 
verging  section  of  the  nozzle.  Contoured  vanes  may  reduce 
the  total  pressure  losses  and  likelihood  of  the  flow  from 
one  end  inter ferring  with  the  flow  from  the  other.  Until 
there  is  m're  experimental  information  {either  cold  flaw  model 
studies  or  motor  firings),  there  is  no  basis  for  expecting 
that  increasing  the  loading  density  of  the  Baseline  Design  by 
decreasing  ADOrt  will  produce  an  overall  increase  in  perfor¬ 
mance. 

A  scries  of  F  vs  t  predictions  were  made  in  which  the 
Baseline  Design  ignition  conditions  were  perturbed.  For 
example,  the  igniter  mass  flow  rate  was  varied  ±  50%  and 
the  igniter  heat  flux  was  varied  ±  25%  while  holding  the  net 
mass  and  heat  generated  by  the  igniter  constant.  The  net 

*  The  lighter  mot^r  case  resulting  £rom  operating  at  a  nominal 
pressure  of  10,000  psi  increases  hy  14%.  However,  the 

maximum  Mach  number  in  the  port  is  .increased  from  0.21  to  0.30, 
which  greatly  increases  the  uncertainties  associated  with 
turning  the  port  flow  90  degrees. 
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- -  BASELINE  DESIGN,  TABLE  5 

- -  SAKE  AS  BASELINE  E25CEPT  1JCKISAL  GPERnil&i 

PRESSURE  IS  REDUCED  FRCJi  15 , OCC  TO  10,000 
?SI  (SEE  PIG.  36) 
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Pig.  38  Calculated  thrust  versa*  tiae  progress  showing  how 
a  baseline  configuration  is  alters  by  reducing 
nominal  operating  pressure. 
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— -  BASELINE  DESIGN,  TABLE  5  _ 

- SAME  AS  BASELINE  EXCEPT  CROSS-SECTIONAL  ARsA  OF 

PORT  DECREASED  BY  50*  {PROPELIANT  WEIGHT  IS  HELD 
CONSTANT) 


Fig,  39  Calculated  thrust  versus  time  programs  showing  how 

a  baseline  configuration  is  altered  by  cross-sectional 
area  of  port  (i.e,,  volumetric  loading  density) 
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-  BASELINE  DESIGN , .  TAKE  5 

_  SAME  AS  BASELINE  EXCEPT  PRESSURE  OF  NOZZLE 

CLOSURE  REMOVAL  IS  DECREASED  FROM  14,000  PSI  TO 
5,000  PSI  AND  CROSS-SECTIONAL  AREA  OF  PORT  IS 
DECREASED  BY  505b. 


Fig.  40  Calculated  thrust  versus  tise  prograss  showing  how 
a  baseline  configuration  is  altered  by  c cabined 
effects  of  reducing  pressure  of  nozzle  closure 
removal  and  reducing  port  area. 
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result  vas  a  variation  in  the  centroid  oi 
than  0.0003  seconds.  These  results  again  demonstrate  that 
the  Baseline  Design  is  operating  in  a  region  where  it  is 
relatively  insensitive  to  off-design  conditions  and  manu¬ 
facturing  variations.  The  strong  incition  sequence  and 
high  nozzle  closure  pressure  (>  1000  psi)  are  key  elements 
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e  transients  of  motor  ignition  processes. 
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under  ignition  can  produce  overpressures  following  nozzle 


closure 

removal.  Such 

a  situation  will  be 

when 

Fig.  43 

is  discussed.) 
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once  a  good  design 
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rrived 
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results  are  xarcely  uninteresting . 

Cue  exams 

le  of 

a  possible  difficult  to  diagnose  experimental  condition  is 
shown  in  Fie.  41.  rioere  41  shows  the  results  of  us inc  an 


igniter  that  does  not  completely  ignite  the  propellant  sur¬ 
face  prior  to  removing  the  nozzle  closure.  Since  the  flame 
spreading  is  not  complete  when  the  closure  blows  out,  the 
charbe r  pressure  can  net  he  sustained  arc  rapidly  decreases 
until  sufficient  burning  surface  is  ignited,  i&iie  the  con¬ 
dition  o?  Fic.  41  is  semewhafc  extreme,  it  is  intended  to 
illustrate  what  sight  happen  if  a  single  weak  igniter  C in  one 
head-end)  is  used  in  a  center  rented  chamber.  An  improperly 
designed  single  igniter  can  have  sufficient  mass  generation 
capacity  but  might  not  have  a  proper  heat  fiur  distribution 
for  rapid  flame  spreading. 


A  series  of  ?  vs  t  calculations  were  made  for  doable 
base  propellants  to  determine  the  influence  {as  compared 
to  the  cccsr-osite  propellant  in  the  Baseline  Design)  of 
higher  burning  rate  exponent  Cn)  and  temper  a  t  ur  e  sensitivity 
of  burning  rate  (c^)  coupled  with  lower  burning  rate  {see 
Table  2? .  These  aBe  the  propellant  changes  that  most  strongly 
aggravate  dynamic  burning  rate  overshoots-  Figure  42  shews 
that  the  higher  burning  rate  overshoots  lave  a  snail  effect 
or  F  vs  t.  Even  though  the  sax  ideas  dynamic  burning  rate 
effect  of  the  double  base  propellant  is  sore  pronounced  and 
sore  sustained  when  compared  to  the  composite  propellant  - 
the  period  of  the  effect  is  essentially  over  before  the 
oczsle  closure  is  removed  {see  Fig-  2S)  and  thus  there  is 
little  influence  on  F  vs  t .  By  constructing  an  off-design 
condition  that  accelerates  the  chamber  pressurization  pro¬ 
cess  and  opens  the  nozzle  before  the  period  of  dynamic  burn¬ 
ing  is  over,  we  can  observe  the  very  detrimental  p  con¬ 
dition  in  Fig.  43.  It  should  be  noted  that  cnlike'*utSe  Baseline 
Design  the  F  vs  t  of  the  Fig.  43  conditions  will  be  greatly 
influenced  by  small  changes  in  parameters  such  as  the  igniter 
sass  flux,  igniter  beat  flux,  A  , /A, .  and  nozzle  closure 
removal  pressure.  Accordingly , if  the  pressure  spike 


*  An  equivalent  interpretation  of  the  cause  of  the  pressure 
overshoot  is  that  the  period  of  boating  prior  to  full 
ignition  ?,s  too  long.  As  a  result,  a  relatively  thick 
preheated  propellant  region  is  burmed  at  an  accelerated 
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-  BASELINE  DESIGN.  TABLE  5 

-  SAME  AS  BASELINE  EXCEPT:  1)  COMPOSITE  PROPELLANT 

IS  REPLACED  BY  DOUBLE  BASE  ^ROPELIANT,  2)  NOMINAL 
OPERATING  PRESSURE  IS  REDUCED  FROM  15,000  TO 


Fig.  43  Calculated  thrust  versus  time  programs  showing  how 
a  baseline  configuration  is  altered  by  a  double  base 
propellant  operating  in  an  of f-aes ign  condition. 
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following  nozzle  closure  removal  is  eliminated  by  small 
adjustments  to  the  motor  and  igniter  parameters,  the  re¬ 
sult  would  be  a  motor  that  is  very  sensitive  to  small 
manufacturing  variations. 

The  results  shown  in  Figs.  35  through  43  illustrate 
that  the  design  (as  it  applies  to  internal  ballistics) 
presented  in  Table  5  represents  a  good  set  of  operating 
points  which  are  relatively  insensitive  to  motor-to-motor 
variations.  It  is  interesting  to  note  that  we  have  used 
analytical  techinques  that  can  predict  pressure  and  burning 
rate  overshoots  to  isolate  operating  domains  where  those 
particular  transient  effects  are  not  detrimental.  Thus 
we  end  up  with  the  ironic  situation  of  describing  a  motor 
whose  performance  can  be  approximately  predicted  without 
considering  dynamic  burning  rate  effects. 


CONCLUSIONS  AND  RECOMMENDATIONS 

Attention  was  focused  on  developing  analytical  techinques  for 
predicting  the  performance  of  impulsive  thrusters  and  on  analyzing 
several  of  the  design  considerations.  Since  specifications  for 
the  Hit  system  are  still  undergoing  modification,  it  was  not  mean¬ 
ingful  to  study  in  depth  designs  to  meet  a  particular  specifica¬ 
tion.  Instead,  typical  designs  were  evaluated  to  indicate  the 
general  trends  and  the  important  system  considerations.  One  of 
the  major  accomclishments  of  the  study  is  the  development  of  ana¬ 
lytical  techniques  which  can  be  readily  applied  to  analyze 
particular  motors. 

The  mathematical  model  for  analyzing  impulsive  thruster  in¬ 
ternal  ballistics  and  performance  contains  several  features  that 
were  previously  unavailable.  In  particular,  the  model  considers 
interact rons  between  chamber  gas  dynamics  and  burning  rate  dynamics 
over  the  entire  duty  cycle  of  the  motor.  Through  the  use  of  the 
Zeldovich-Novozhilov  method  transient  burning  rates  can  be  calcu¬ 
lated  from  input  values  that  are  more  readily  obtainable  than  the 
input  values  required  by  previously  reported  models. 

One  of  the  explicit  results  of  this  study  is  the  impulsive 
thruster  design  of  Table  5.  The  overall  quality  and  practicality 
of  the  design  can  be  evaluated  from  two  aspects:  1)  the  motor 
hardware  and  2)  the  internal  ballistics.  From  the  aspect  of 
motor  hardware  the  design  may  be  idealistic  since  it  calls  for  the 
fabrication  of  high  strength  case  joints,  nozzle  to  case  flanges 
and  igniter  bosses  that  use  difficult  to  work  with  steels.  Accord¬ 
ingly,  W  ;/W  requirement  in  the  specification  (Table  3)  will 
be  difficult  to  meet.  From  the  aspect  cf  internal  ballistics, 
the  design  looks  very  good.  The  design  was  placed  in  an  operating 
domain  where  generous  manufacturing  tolerances  and  off-design 
conditions  car  be  permitted  without  serious  performance  losses. 

The  use  of  a  strong  igniter  and  nozzle  closure  designed  to  open 
above  1000  psi  are  key  items  in  obtaining  control  over  the  ignition 
interval  and  a  reproducible  centroid  of  the  thrust  versus  time 
program. . 

As  a  result  of  this  study,  we  have  formed  several  opinions 
ccr.c  err.inc  how  the  propulsion  system  development  should  procet.d.  We 
believe  that  canter-vented  impulsive  thrusters  can  be  built  with 
good  oer  formar.ee .  However,  proceeding  with  the  design  of  a  motor 
to  meet  the  requirements  of  Table  3  without  concurrent  experimental 
evaluations  may  delay  obtaining  information  needed  to  size  the  moto 
components,  (e.g.,  thrust  vector  alignment  of  center- vented  motors, 
effect  of  low  A  on  repeatability,  and  plume  interference). 

The  dancer  is  tRStuthi  propulsion  problems  (and  subsequent  fixes 
using  heavier  motors  to  achieve  the  desired  a$)  will  be  discovered 
sequentially  during  future  development  programs.  We  would  like  to 
focus  attention  on  items  that  may  require  a  quantitative  determina¬ 
tion  to  upgrade  the  propulsion  systems  that  are  being  developed. 
7hes-->  ite~s  are: 

I ■  what  is  a  realistic  upper  limit  of  propellant  loading 
:crs:ty  in  center-vented  motors  without  destroying  the 
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desired  characteristics  hy  incieasmq  t crust  nis- 

alignmentkanc!v',fail-of  f . 

2)  What  is  the  effect  of  out-gassing  from  plastic  cases 
(if  used)  and  case  insulation. 

3)  To  what  extent  do  impinging  internal  flows  cause  inter¬ 
actions  that  produce  nonrepeatabie  and  misaligned 
thrusts . 

4)  To  what  extent  can  back  flew  of  the  plume  be  alleviated 
by  design  changes  involving  such  items  as  propellants, 
pressure  level,  nozzle  contouring  (cross-sectional  and 
axial  contours),  and  exit  cone  surface  coatings. 

The  effects,  characterization,  and  alternatives  resulting  from 
the  above  items  are  included  in  summary  in  Table  7. 

The  magnitude  of  the  effects  associated  with  Items  1 
through  4  can  be  evaluated  only  by  direct  testing.  Kote  that  the 
tests  do  not  require  that  the  ter.-:  motc_s  use  elaborate  fabrica¬ 
tion  techniques,  high  safety  factor  hardware  will  do.  Analytical 
treatments  associated  with  each  item  would  be  valuable  in  finding 
the  best  operating  points  and  design  fixes  (if  required) .  In 
particular,  we  strongly  recommend  studies  be  conducted  on  the 
following  internal  ballistics  and  gas  dynamics  items: 

1)  Diagnostic  tests  under  vacuum  conditions  to  evaluate 
those  designs  and  propellants  that  alleviate  back-flow 
of  the  plume  (and  not  merely  measure  the  plume  spillage 
for  a  few  configurations) . 

2)  Systematic  tests  in  a  high  precision  multi -component 
thrust  stand  (or  the  equivalent  ballistics  device) 

to  quantitatively  establish  the  influence  of  internal 
flows  on  the  thrust  vector  alignment. 

3)  Conduct  a  test  series  to  establish  the  repeatability  of 
the  thrust  versus  time  trace. 

Finally,  we  recommend  that  the  analytical  techniques 
developed  during  this  study  be  subjected  to  the  test  of  predict¬ 
ing  the  performance  of  oarticular  developmental  motors  for  the 
Hit  systems. 
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